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Introduction 


Asymptotic Modal Analysis (AMA) is a method which has been used successfully 
to model and understand linear dynamical systems with many participating modes. The 
AMA method was originally developed by Dowell [1] to show the relationship between 
statistical energy analysis (SEA) and classical modal analysis (CMA). In the limit of a large 
number of modes of a vibrating system, the classical modal analysis result can be shown to 
be equivalent to the statistical energy analysis result. As the CMA result evolves into the 
SEA result, a number of systematic assumptions are made. Most of these assumptions are 
based upon the supposition that the number of modes approaches infinity. It is for this 
reason that the term "asymptotic" is used. AMA is the asymptotic result of taking the limit 
of CMA as the number of modes approaches infinity. AMA refers to any of the 
intermediate results between CMA and SEA, as well as the SEA result which is derived 
from CMA. 

The main advantage of the AMA method is that individual modal characteristics are 
not required in the model or computations. By contrast, CMA requires that each modal 
parameter be evaluated at each frequency. In the latter, contributions from each mode are 
computed and the final answer is obtained by summing over all the modes in the particular 
band of interest. AMA evaluates modal parameters only at their center frequency, and does 
not sum the individual contributions from each mode in order to obtain a final result. The 
method is similar to SEA in this respect. However, SEA is only capable of obtaining 
spatial averages or means, as it is a statistical method. Since AMA is systematically derived 
from CMA, it can obtain local spatial information as well. 


Summary of Research Conducted Under the Grant 

AMA is considered in this work within the framework of those linear dynamical 
systems called structural-acoustics. Some general analytical results were obtained for the 
case where part of an otherwise rigid enclosure is vibrating with a large number of 
structural modes, which excite either a finite or an infinite number of acoustic modes in the 
interior. In Reference 2, it was shown that Statistical Energy Analysis is an asymptotic 
limit of classical modal analysis. The basic asymptotic theory for structural wall-acoustic 
cavity interaction was described, and several numerical examples were presented for the 
acoustic cavity response. 

In Reference 3, the AMA results for an infinite number of structural modes were 
applied to the specific case of a rectangular acoustic cavity. The rectangular acoustic cavity 
had rigid walls except for a flexible portion on one of its walls. The flexible portion 
vibrated with an infinite number of structural m od es which were temporally uncorrelated 
("white noise"). The interior sound field was studied using asymptotic modal analysis. 
Both local mean-square pressures and spatially averaged mean-square pressures in a 
bandwidth were computed. The AMA results were compared to CMA results. In addition, 
the effect of varying the size and the location of the flexible vibrating portion was studied. 

It was found that the interior sound pressure levels are nearly uniform, with 
exceptions occurring at the boundaries. The boundaries exhibited elevated levels which 
were 8, 4, and 2 times greater than the interior level, for comers, edges and walls 
respectively. These areas of elevated sound pressure level were termed "intensification 
zones." 

In the numerical study, the size and location of the flexible portion of the wall was 
varied, and the resulting effects on the interior sound field were studied. An interesting 
phenomenon was discovered. When a small source (smaller than an acoustic wavelength) 
was placed on a point of symmetry of the cavity, elevated sound pressure levels occurred 
along the planes of symmetry, creating new effective boundaries of the sub-divided cavity. 


These effective boundaries retained the same ratios of 8, 4, and 2 for differences between 
comer, edge, and wall intensification, relative to the interior (of each sub-cavity). Another 
result of the numerical study was that sound sources placed in the comer, rather than 
elsewhere on the wall away from a boundary, raise the sound pressure levels by a factor of 
four. It was further deduced that sound sources placed on an edge would raise sound 
pressure levels by a factor of two. 

In Reference 4, the intensification zones were studied using AMA and other 
traditional methods. The system parameters which define the structure of an intensification 
zone were determined. Using AMA and non-dimensionalizing the spatial variable by the 
center frequency wavenumber, the intensification zone description was found to be 
independent of the cavity dimensions. In fact, the non-dimensional spatial variation of 
sound pressure level was only dependent upon the ratio of bandwidth to center frequency 
(fb/fc)- Upon further study, it was shown that the bandwidth effect is second order in 

fb/fo 

The independence of cavity dimensions suggested that intensification could be 
modelled as a local system. This led to modelling using oblique incidence sound waves. 
This alternative method gave identical results. Reflections from an absorptive surface were 
also studied and an analytical expression for the sound pressure levels near an absorptive 
wall was developed. This result can be expressed as the sum of a rigid wall term plus a 
term containing the impedance information for the absorptive surface. 

Experiments were performed in order to validate the AMA method for the 
structural-acoustic application. The results of the experiments are presented in Reference 5. 
The experiments tested for uniformity of the interior, intensification near comers, edges 
and walls of the cavity, and accuracy of the method for predicting sound pressure levels in 
the interior. In all cases, the agreement between experiment and analysis was very good. 

Qualitatively, it was shown from the experimental results that AMA works best for 
moderate bandwidths. If the bandwidth is too narrow, there may not be enough modes for 



the AMA assumptions to apply. However, if the bandwidth is too wide, there will be 
enough modes, but there may be too much variation in the modal variables which are 
assumed constant throughout the band in the AMA methodology. 

The AMA method was used to predict the spatially averaged sound pressure levels 
in the cavity interior from the measured acceleration of the vibrating wall. The predicted 
levels were compared to the experimentally measured levels and good agreement was 
obtained. 

Further theoretical work was done to develop an AMA result which is valid for the 
case where the number of structural modes is finite and the number of acoustic modes is 
infinite. These results are reported in Reference 6. The previous AMA work considered 
the case where the number of structural modes was very large or, formally, infinite . It 
was found that different AMA results are obtained if the limits of an infinite number of 
structural modes or an infinite number of acoustic modes are taken in a different order. 
Therefore, when there are both a large number of acoustic modes and a large number of 
structural modes, it is important to use the appropriate limiting sequence. This decision 
should be based on which set of modes is approaching infinity most rapidly. 

References 2 through 6 are presented in Appendices A through E of this document. 

Recommendations for Future Work 

One outstanding issue is worthy of mention. In the usual linear dynamical systems 
theory, be it CMA, SEA or AMA, the assumption of small damping is often made. From 
the experimental work done under the present grant, there is clear evidence that all of these 
theoretical models need to be extended to describe accurately the physical phenomena for 
large damping. The high damping case, it may be noted, is of considerable practical 
importance. 


References 


1 . Dowell, E. H., "Vibration Induced Noise in Aircraft: Asymptotic Modal Analysis and 
Statistical Energy Analysis of Dynamical Systems," R82-1 12447, United 
Technologies Research Center, 1983. 

2. Kubota, Y., H. D. Dionne and E. H. Dowell, "Asymptotic Modal 

Analysis and Statistical Energy Analysis of an Acoustic Cavity," Journal of Vibration. 
Acoustics. Stress and Reliability in Design . 110 , 371 - 376, 1988. 


3. Peretti, L. F. and E. H. Dowell, "Asymptotic Modal Analysis Of A Rectangular 
Acoustic Cavity Excited by Wall Vibration," AIAA Journal. 3 0. 1 191-1198, 1992. 

4. Peretti, L.F. and E.H. Dowell, "A Study of Intensification Zones in a Rectangular 
Acoustic Cavity," AIAA Journal . 30, 1 199-1206, 1992. 

5. Peretti, L. F. and E. H. Dowell, "Experimental Verification of the Asymptotic Modal 
Analysis Method As Applied To A Rectangular Acoustic Cavity Excited By Structural 
Vibration," presented at the ASME Winter Annual Meeting, November 1990 in 
Dallas, TX, to appear in ASME Journal of Vibrations and Acoustics , 1992. 

6. Peretti, L. F., "Chapter 4: Asymptotic Modal Analysis for Structural-Acoustic 
Systems with Finite Structural and Infinite Acoustic Modes," Asymptotic Modal 
Analysis for Structural-Acoustic Systems . Ph.D. Thesis, Duke University, 1991. 


Appendix A: 


Kubota, Y., H. D. Dionne and E. H. Dowell, "Asymptotic Modal Analysis and 
Statistical Energy Analysis of an Acoustic Cavity," Journal of Vibration. Acoustics. Stress 
and Reliability in Design . 110 , 371 - 376, 1988. 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Y. Kubota 2 

Visiting Scholar. 

H. D. Dionne 

Research Assistant. 

E. H. Dowell 

Dean. 

cnt of Mechanical Engineering 
and Materials Science, 
School of Engineering, 
' Duke University, 

Durham, North Carolina 27706 


Asymptotic Modal Analysis and 
Statistical Energy Analysis of an 
Acoustic Cavity 1 

One of the outstanding theoretical questions in interior noise is the connection be- 
tween modal analysis and statistical energy analysis. Recently substantial progress 
has been made in understanding this connection for structural vibrations including 
both fundamental theoretical work and experimental verification . It has been shown 
that many of the results of Statistical Energy Analysis can be derived as an asymp- 
totic limit of classical modal analysis and thus this approach is called Asymptotic 
Modal Analysis. The basic asymptotic theory for structural wall-acoustic cavity in- 
teraction is described in this paper. Several numerical examples are presented for 
acoustic cavity response with a prescribed wall motion to illustrate the key results of 
the theory . 


faction 

[1] showed that the results commonly referred to as 
al Energy Analysis (SEA) [2] as developed by Lyon 
can be obtained by studying the asymptotic 
of classical modal analysis (CM A) for a general, 
[structural system; those asymptotic results are called 
: Modal Analysis (AMA). In [1], moreover, specific 
ations were made for structural-acoustic systems and 
ling subsystems. Since AMA results can be derived 
atically from CMA, AMA allows an assessment of the 
ptions and consequent simplifications which are made 
ain such results and a combination of CMA and AMA 
K) may prove useful in applications. 

3], a comparison of AMA (or SEA) and CMA was made 
response of a single general linear structure and the 
Iptotic characteristics of AMA were discussed. It was 
fjp that the asymptotic behavior of AMA (or SEA) 
s upon the number of modes in a frequency interval of 
iltt and the location of the point forces. Moreover, 
^toticaily, all points on the structure except for some 
« ones have the same response; the exceptional points are 
F the points of excitation and near the structural system 
*&ry. Some numerical examples for a beam were 
Ated in [3]. 

'Hi the response of a rectangular plate under a point ran- 
1 force was investigated experimentally and theoretically, 
^were two objectives in this study. The first one was to 
^^trate experimentally the manner in which the asymp- 
: Emit is approached. The second one was to show ex- 
•entally that the response of almost all points of the plate 
the same in the asymptotic limit. 

^rimental measurements were carried out for several 

w °rk was supported, in part, by NASA Research Gram NAG- 1 -709 
P* Structural Acoustics Branch of the Langley Research Center Dr Kevin 
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or Vibration, Arorvnrs, Strfss, ani> RnuBimv in Dimov 
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frequency bandwldths in which the maximum one corresponds 
to 47.1 modes of plate vibration, for several band center fre- 
quencies, and for two locations of the point force which were 
the center of the plate and near the plate edge, respectively. All 
experimental results verified the asymptotic behavior of the 
plate response which is predicted by AMA and were in good 
quantitative agreement with the theory as the number of 
modes became sufficiently large. 

The important work of Crandall and his colleagues should 
also be cited [5-8]. Their work has been directed primarily 
toward those aspects of structural response under wide-band 
random excitation where SEA-type results of the usual sort do 
not hold. Both analytical and experimental work has been car- 
ried out. 

Modal sum and image sum analytical methods have been 
used by Crandall, The modal sum methods are similar to those 
employed here. Image sum methods involve repeatedly using 
traveling waves on an infinite structure to form an appropriate 
cancellation of image sources to model finite plate response. 
Both methods ha ve led to the identification of “intensification 
zones” where the local response is higher than the nearly 
uniform response of the rest of the structure. Special attention 
has been given to square and circular plates where such zones 
are most pronounced. 

Here the emphasis is on extending the earlier work in the 
area of structuraFacoustic systems. A theoretical discussion of 
this subject is given next, followed by a section with numerical 
examples. A complementary discussion of this subject for 
acoustic sources internal to the acoustic cavity is contained in 
[9). 

Theoretical Discussion: Noise Transmission Through a 
Structural Wall into an Interior Acoustic Cavity in the 
Asymptotic Modal Analysis (AMA) Limit 

The structural wall modal equation of motion is (see [1], |3], 
and [4] for technical background and nomenclature). 
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M m \Q m + Qitt 

where the modal expansion Tor the wall deflection is 

m 

and the structural generalized mass is 


(1) 


(2) 


M, 




mA m 2 dxdy 


and the generalized force due to a given external pressure is 

Qm = jj P E ^ m dxdy (3) 

The acoustic cavity modal equation is 

P, + 2t*a?P r +af 2 p,=Q? w 

where the modal expansion for the acoustic cavity pressure is 

P r (t)F,(x,y,z ) 

P = Po4L JfA 


ORIGINAL : 'S 
OF POOR i'V 


This is the basic expression for the power spectra of the cav** 
pressure in terms of the power spectra of the wall acceleration 
I. When the number of excited structural modes is 1*^ 
AM — oo, it can be shown that 

This means that the power spectra of the wall respond * 
uncorrelated in space. This assumption is reasonable for 1** 
AM , because 

1 


0 


(X*x\ y*y* 


AM 


<Mu >\x,y y x\y*)- 


UJ* 


(5) 


and the acoustic generalized mass is 

M A ■ -y j j j Fjdxdydz (6) 

and the generalized acceleration due to the structural wall is 

Define /, a nondimensional cavity pressure. 


f(t,x,y,z)> 


PqCq 


constant (x = x\ y^y*) 
as AM— oo 

Recall [1,3] 

Yd u>) 

m n 

u 2 H n (-u)- £ £ ^ m (x i ,y,)i n (x’ J ,y , J )* Fij {u>) 

f j 

(13) is readily derived from the above relationship and invofc- 
ing the basic methods of AMA. $ F is the cross spectra lot 
point forces, F n F Jt on the structural wall. Also foe a 
smoothly varying power spectrum, it is assumed that 

**(tt)s* # (« f ) 

This is just the usual white noise assumption. Thus, equaiiot 
(11) becomes 

a f - , v v' v' F r ( x >y,z) 

^/(w;x^,z)=-p— 4>w(^c) L L — JjA — 


From (5) the auto-power spectrum of/ may be determined as 
« , v F r (x,y,z) F s {x,y,z) f , ^ 

tyu.JJ'.z)- 2j Zj jjt Jjz *P,P S ( 7 ) 

where the cross-spectra are defined as 

^ p r p s m ~\^ a} R p r p s ^ eiul7dr 

and the cross-correlations of the modal generalized pressure 
coordinates are 

R Pr P ~ Hm A;\\pAt)P,V + T)dt 
Similarly from (6) the cross-power spectra of Q? and Q* are 

V«r m -pr\ L S L, 

( o) \x,y >x*j* )dxdydx* dy * (8) 

From (4) and standard random response theory, the rela- 
tionship between <b P P and $ is 

r S Vr Vs 

$ PrPs (u) -H A ( -V^qWqW (u>) (9) 

where the modal transfer function is defined as 

H? <u) j — — (10) 

r-w 2 +as f +2 itfufu] 

From (7), (8), and (9) 

1 r-i F r (x,y,z) F s (x,y,z) 

* / (w;x,y,z)=— r Ij L—JjA ma 

F,(x,y,z 0 )F s (x*,y*,Zi 0 *W„ 
(uSAr ,y,x’,y*)dxdydx*dy* (11) 


V 1 

F, (x,y,z) 


.//? («)•//? (-«) 


• | F, (x,y,z 0 )F s (x,y,z 0 )dxdy 


(IS 


The mean square response of the nondimensional cavfej 
pressure is 

^ I 

M r uy ^ 

Note that (16) is the result for a large number of strurlrfH| 
modes, AM — oo, and a moderate number of acoustic c 
modes, AN a . i jf : T- 

Taking a spatial average of (16), and noting that M, 
fj? , <Fj> do not vary rapidly with respect to modal nutaW^ 
r, for large AN a , equation (16) becomes 


<p l > 

(Po^o) 


T A 

T 


F 1 r (x,y,z 0 )dxdy 


*F 


where < > denotes spatial average. 

Now r consider a cavity acoustic modal function 

F,{x t y,z)=X r (x)Y,(y)Z,(z) 

Take the plane at z = z 0 as the boundary of the a ? 0US **jL 
where the structural wall is vibrating. Z,(z 0 ) is 
dependent of mode number r or it can be so normalized, 
for large AN A , 
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. , - Vj ^ ^ V 2 ^ 7:/^ \ Equation (25) will be dominated by the terms for which r = s, 

F*(xj,z 0 )dxdy= ^A F <K> a f < = j Thus equa(ion (2 / } becomes 


-~A f »AN a 


<Z?> 

</7> 


(19) */(<<W-,2)s-i- £ £ E 


•w 4 l// m (w) i 2 «<Jy (w) 


M A = <F;> 

<*' 2 > Aw“Aw $*(w c ) 




gp.- j r r 

K> denotes a spatial average. <F C 2 > is a volume Am E J 

of A; 2 and <F c 2 > a is an area average over the 

E structural wall. Hence, equation (17) becomes as mean square response for the nondimensional cavity 

acoustic pressure is 

JLjL-^_(^V <ga ~ ,20) 

tWF" < u y !t<z;> J " 


tfom the AMA results for structural wall motion, 

<^5>. ~ ^ 

I **“ 4 Au) JW*«& 

then, equation (20) becomes 
p> / r \ 2 AM AN A / A F \ 2 Aw 4 

Ln) \ 4 / Aw Aw' 1 \ K / Aw 


75 r* . 

? *T^F" a )o */<-***>** 


Using (26), 

(2i) •^ £- 7 __ prE E E 


&rm 

m i u y m U A * f m 


(22) where 


w c <F>Aw 


a rm m fm w m + fr* (£*m w m + ^ ) 

^rm a (^m“W^ 2 ) 2 ^4w m ^(f m CJ m + ^^)(f m ^ + tfw m ) 

This result assumes that the acoustic and structural resonances 
are well separated. 

(23) See Crandall and Mark [10], p. 72, for the evaluation of 


/» * 

f ie result for AM — oo and AN a — oo. (16) is the result ] 0 ^ A (w)\ 1 \ H m (w)\ 1 dw 

- oo and finite AN a . 

bove derivation is not the only one which leads to (23). When fm* are small, then 

instructive to consider an alternative derivation which a rm s f m wj, + f^w? 3 (28) 

ben equation (12) can not be assumed, then When Aw/w c , Aw /u f << 1, then 

' V V F,(W> F,0w> »„ + «>■ (29) 

~ ^ r 5 A/? * A/^ (<*>) All (29) really assumes is that w m — is small. Hence in sum- 

ming, w m and <u ) A , one need not distinguish between them, 
r (• ci* Taking a spatial average, and noting that 

W), lW ) L/^MoWAx-'S.Zo-)** (24) 

[ do not vary rapidly with modal number for large AM and AN a 

r (u;x,y,x m t y*)dxdydx*dy* (II) and letting 

w f = w A , Aw = Aw' 4 , 

Utuvj V T f / w / • 2 ej / t then equation (27) becomes 

K^ ) .SLww»,wv*.(.i <?> . , ^ 

L . . ^ ^ ‘ 4 <" T ' 


( w ;*,>> , x 1 * ,y * )dxdydx *dy* 


m n 

3 Hj - w )* E E 


*«utc equation (24) into equation (11), then 

« X . )= _L T V V V* F r^.z) F, (xp\z) where 

\n Lj Lj Li Ls * fcjA 

r s m n iVI r JVI s 


(«r) E E 

r m Prm 


(u)H A ( — u)»o) 2 H m (w)w 2 //„ ( — w) 

•E E 


IL F r(x,y,z 0 )F s (x , % y*,z , 0 


)'l‘ m (x,yW„ 


( x* t y * )dxdydx * dy * 


°r *^(fr + fr ) 

For large AM , AA^ . we can write 

E E - 3 — E E ( 3 i) 

r m Pfm Pf r m 


where l/0 t is the average of (l/0 rm ) 

r m 

(25) \/0 c can be approximated by an integral expression 
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2' x 3' x 7 ' Cavity 

2' x 3 1 Flexible Vtbrating Wall 

Fig. 1 Cavity geometry 
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Thus, equation (30) becomes 

<p 2 > _ ( * \ 2 1 ^ <Fj > yy 

(Pod,) 2 ~ ' 4 / Au) F 5 MlMf 1 01* [ ( if - F " “ f) 

• E I] d m {X„y,)dm 

r m 

Now, the summation ^ ^ can be written as 

r m 

Srm S T, E 

r m 

= E £ ^(W.)[JL F d x >y^n,(x,y)dxdy\ 

r m F 

* [11^ ^ (jf ‘ u 1 * ’ z ^ m (x ‘ u'*) d **4>'*] 

,y\zo) 

•[ E i 2 m(x,,y l H m {x,y)^ m {x , o'’)]dxdydx t dy' (Mj 

L m 

For large AM, the summation* ^ ^ (at, v) 




(x*, y*), will be dominated by x = x\ y = y* 
and 


« K»Tt 
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Fig. 2 
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as AM— oo 

Thus, S rm can be written 

S rm a AMI < > l 2 • E j ^ P, (x,y,z a )dxdydx'df 

= AMI < > ] 2 'A f E ( ( ^ (x,y,z 0 )dxdy 

r J ^ a F 

For large AN 4 , S rm becomes 


(Jii 


oe 


S rm a AM-AN a A 1 f [ <tf>] 2 ' 


<n> 


7 ' « 3 i 7 ' CAVITY 


<Z?> 

Substitute equation (37) into equation (33) and note that 

M r = M P <*p 2 > 

M 4 = <F^> 


HUMBER of HOOES vs CENTER FREQUENCY 

Fig. 3 Number of modes versus center frequency 
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Then equation (33) becomes 
<p 2 > ( t \ 2 AM AN A / A F \ 2 


(32) 


I I f“ 1 
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Mp u c£ctc <Z}> 

Equation (38) is the same as equation (23) for w f = w^ 4 . ^ ^ 
Aw* 1 . 
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Illustrative Numerical Examples 

: Consider a rectangular cavity of dimensions 2 ' x 3' x 7\ 

r** Fig. 1. The 2' x 3' wall is taken to vibrate in a large 
pwnber of structural modes, AM — oo.We wish to compare 
U* results of Asymptotic Modal Analysis for the acoustic 
Ovity response, equation (20), to that of Classical Modal 
Analysis, equation (16). For simplicity and to focus the discus- 

r? n on acoustic response, the wall motion is assumed to be 
liven. 

^The ratio of the spatial average of equation (16), CM A, to 
•Oration (20), AMA, is shown in Fig. 2. This ratio is shown 
ET three different frequency bandwidths and as a function of 
H* center frequency. The bandwidth, Aa>, is defined as 

[ 

the center frequency, w c , as 

w min are the maximum and minimum frequencies 
the frequency interval. All acoustic modes are assumed to 
v * the same modal critical damping ratio, , for simplicity, 
can be seen from Fig. 2, all results approach unity as the 
ler frequency becomes large. The approach is modestly 
rc rapid for the smaller bandwidths. To better understand 
results of Fig. 2, consider Fig. 3. Here is shown the number 
°des, AM , in a given bandwidth as a function of center 
Wncy. It is expected that CMA/AMA — 1 as AM — ce. 


Fig. 7 CMA/AMA ratio versus center frequency 

This clearly explains why CMA/AMA - I as the center fre- 
quency becomes large. However considering Fig. 3 it is 
perhaps initially surprising that in Fig. 2 the asymptotic ap- 
proach is modestly more rapid for the smaller bandwidths. 
This is explained by recalling that in AMA we neglect the 
variation of natural frequencies from one mode to another 
within a given bandwidth. Thus while increasing the band- 
width increases the number of modes (which improves the ac- 
curacy of AMA), it also increases the variation among the 
natural frequencies (thus decreasing the accuracy of AMA). 
The net effect of these two competing factors is to decrease 
modestly the accuracy of AMA as the bandwidth is increased. 

The above results are all for the spatially averaged pressure 
inside the cavity. Now consider the ratio, CMA/AMA, for 
local response at various positions within the acoustic cavity. 
In Figs. 4, 5, 6 the ratio is shown versus position in three dif- 
ferent diagonal planes in the cavity. A single frequency inter- 
val, 1400-1800 Hz, is considered. As can be seen the ratio is 
everywhere nearly constant and close to unity, except along a 
cavity side, edge and, especially, at a cavity corner. On 
theoretical grounds the ratio is expected to be, asymptotically, 

side: 2 

edge: 4 

corner: 8 

To see how these expected asymptotes are approached, con- 
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sider Fig. 7. One corner point, A, and three edge points, B, C, 
D, arc considered. The ratio, CM A/ AM A, is shown as a func- 
tion of center frequency for a fixed bandwidth. As expected 
the ratio approaches 8 for the corner point and 4 for the edge 
points. 

Conclusions 

Several points were brought out in the text that bear em- 
phasis. When spatially averaged, the Classical Modal Analysis 
(CMA) response approaches the Asymptotic Modal Analysis 
(AMA) response more rapidly as the number of modes in- 
creases, A N a — oo. However, information about local 
response intensification is lost in the averaging process. Ex- 
amination of local response as A N A — - oo revealed an almost 
uniform response throughout the interior cavity field with the 
exception of peaks at the boundaries, i.e., sides, edges, and 
corners. 

Although a larger bandwidth at a given center frequency 
contains more excited modes than a smaller bandwidth, AMA 
is often modestly more accurate in the smaller bandwidth. 
This is because the increased variation of natural frequencies 
from one mode to the next in the larger bandwidth decreases 
the accuracy of the AMA result. Finally it. is noted that all 
AMA asymptotes were approached from below by CMA with 
fixed bandwidth and increasing center frequency. 

Future Work 

• Basic Theory for Asymptotic Modal Analysis 

AA/* number of structural modes 
AN a « number of cavity acoustic modes 

The theory for 

AA/^oo, AN a finite 
and AA/ — oo, AN a — oo 

has been worked out as has been discussed in the text. The 
case AA/ finite and A N A — oo has yet to be considered. Of 


course the case, AA/ finite and A N A finite, is 
classical modal analysis. 

Numerical Studies. Further work might include stud*, 
the reduction of the moving portion of the wall to 
source. In this paper an entire wall is taken as the 
area. It would also be interesting to consider more 
plex cavity geometries. The rectangular cavity shapt^T 
chosen for this study because its modal char act erisiiaTl! 
well known and therefore comparison of CMA 
AMA was a readily feasible task. 

Experimental Studies. An experimental study should W 
undertaken to assess the predictions of the numeric 
studies. A rectangular box with a portion of one 
vibrating might be considered. 
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the location along the wall and the size of the J*' {rom clslssica | modal analysis to that obtained from 

Interior were computed as a ratio of the result ceneral this ratio converged to 1.0 as the 

asymptotic modal analysis for the various plate conr.guratmns^I^^^l^^ ^ ^ (he „ cita ,ion location 
number of responding modes increased lntcns.rication ^ efficient and accurate in solving the 

■■ — - - - ■” -*• di ”” ,c ‘ 

problems with a large number of responding modes. 


Nomenclature 

A = area 
c = speed of sound 

f = cavity acoustic modal function 

/, = cavity dimension 

M = generalized mass 

p = pressure 

r = modal index 

y = volume 

A ,Y,Z = acousdc^odal function component dependent 

on x,y,z 

spatial position coordinates 
number of acoustic modes 
damping ratio 
density 

power spectrum 
frequency 

spatially averaged quantity 
= time derivative 
= rms 

Subscripts 

c = center frequency 

/ « flexible 

o - reference value 

r = acoustic modal index 

w *= pertaining to the flexible wall 


x,y,z 
A AT 


) = 


Superscripts 


acoustic 


Introduction 


Hiuvuwv..*- 

<p\ OUPLED structural-acoustic systems are * t 

Ls often in aeronautical applications. Accurate, ^ 

means of analysis are central to t e esl ®. analytical meth- 
the desired sound transmittal properties. Two 
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_ j, tha , are commonly used to solve such coupled systems are 
classical moS analysis (CMA) and s.a.is.ical energy analysis 
Sfai Recently Dowell 1 and others 2 ' 4 have developed an 
additional method, asymptotic modal analysis (AMA), which 
ran also be applied to structural-acoustic systems. 

AMA incorporates the advantages of the other two meth- 
ods Providing there are a large number of modes, the CMA 
result and the 8 AMA result are nearly identical. Since it does 
not need to take the individual modal contributions into ac 
coum^the^omputation cost of AMA is significantly An 
added advantage of AMA is that the degree of generality in the 
final result "an be controlled by adjusting the types of assump- 
tS and/or simplifications made in the derivation This fu- 
ture allows the use of AMA to obtain results identical to SEA 
or for example, to relax the averaging simplifications and 
obtain local response results of which SEA is not capa e. 

To explore the capabilities of AMA, the interior sound field 
of a rectangular acoustic cavity was analyzed mathematically. 
The rario of response predicted by CMA to that predicted by 
AMA was calculated either as a spatial average or at partwutar 
locations 5 inside the cavity. Five of the cavity walls were rigid 
A random white noise sound 

the sixth wall into the interior of the cavity. The fiewble vibrat 
ine portion was varied in size and location, and the ^suiting 
sound pressure levels (mean square pressures) in the interior 
SteteSEd “ s 'ng AMA and CMA and .ten 

Local response peaks, or intensification zones, were oo 
served at boundary points, while the response m the mterio 
eSon wawiearty uniform. Sound pressure levels were af- 
footed by^ic location and ate oft* 

of the wall (sound source). In general, sound pressure levels 
increased by a factor of 4 when a point was placed m the 
corneT Also for small sound sources (smaller than an acoustic 
wavelength) that are placed in the 

symmetry become regions of increased sound pressure level. 

Background 

SEA has been used to study the high-frequency interaction 
b.»«n large, oompfox, multimodal rate. 
snaces The basic assumption underlying SEA is that tne y 
namic parameters in the system behave stochastically. SEA 

SupTed^^ 
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The main advantage of SEA is its ability to describe the sound 
field without having to consider the individual modes. Statisti- 
cal energy analysis also allows for a much simpler description 
of the system, requiring only parameters such as modal den- 
sity, average modal damping, and certain averages of modal 
impedance to sound sources. The most significant disadvan- 
tage of using a statistical approach is that it is only valid for 
systems whose order is sufficiently high that the stochastic 
assumptions apply. Certain frequency bandwidths may not 
contain enough modes to allow the underlying assumptions to 
hold, rendering the SEA result unreliable. In addition, the 
local response information is lost in the SEA treatment. The 
text by Lyon 5 is the standard reference on SEA. 

Dowell 1 has shown that results identical to those calculated 
using SEA can be obtained by studying the asymptotic behav- 
ior of CMA for a general, linear (structural) system; this 
asymptotic approach is AM A. AM A is basically a modal sum 
method. It possesses all of the computational advantages of 
SEA, in that the individual modal characteristics do not play 
a role in the asymptotic analysis. Additionally, AMA has ad- 
vantages that SEA does not. Since AMA results can be derived 
systematically from CMA, AMA allows an assessment of the 
assumptions and consequent simplifications that are made to 
obtain such results. Also, by using a combination of CMA and 
AMA, results can be obtained for all frequency bandwidths of 
interest, not just those with a sufficiently high number of 
modes. And finally, AMA has predicted local response peaks, 
or intensification zones, results unobtainable using SEA. 3 ’ 4 

Previous work for structural systems has shown that the 
asymptotic behavior of AMA depends on the number of 
modes in a frequency interval of interest and the location of 
point forces. In the limit of an infinite number of modes, all 
points on the structure have the same response except for some 
special areas. The exceptional areas (intensification zones) are 
near the points of excitation and near the structural system 
boundary. 3 - 4 Numerical examples were presented for a beam in 
Ref. 2. Crandall, 6 Itao and Crandall, 7 and Crandall and Kul- 
vets 8 experimentally found intensification zones in their work 
with structures. The response of a rectangular plate under a 
point random force was investigated by Kubota and Dowell, 3 
and AMA calculations were found to agree closely with exper- 
imental measurements. 

Work has also been done using AMA for structural-acoustic 
systems. Kubota et al. 4 examined a rectangular acoustic cavity 
with one vibrating wall (the other five rigid). They assumed 
that the vibrating wall had an infinite number of structural 
modes responding and that the entire wall was oscillating. The 
results obtained from the numerical study indicated that the 
spatially averaged CMA response approaches the AMA re- 
sponse as the number of modes increases. The local asymptotic 
response revealed an almost uniform distribution in the cavity 
interior, with peaks at the boundaries (sides, edges, and cor- 
ners) of the cavity. Here, an AMA method is presented and 
new results obtained for structural-acoustic systems where 
only a portion of the wall vibrates rather than the entire wall, 
and both the size and location of the oscillating portion are 
varied. 


Theory 

Most coupled structural acoustic problems are modeled us- 
ing either CMA, summing for the response of each mode, or 
SEA, which combines the predicted energies of the subsystems 
and coupling loss factors to obtain a final result. In this work, 
a comparison is made between the CMA result and the AMA 
result as the number of acoustic modes and the number of 
structural modes approaches infinity. Note that the spatially 
averaged AMA result is identical to the SEA result. 

In order to calculate the response of the interior acoustic 
cavity to the transmission of noise through a structural wall on 
its boundary, both the structural modes of the wall and the 


acoustic modes of the interior must be considered. As derived 
in Ref. 4, the CMA result for the mean square pressure at a 
point (*,>>, z) in a structural-acoustic enclosure is 


p a i dr* r ^ v F r( x >y> z ) 

( Pocl ) 2 “ 4 V * 


[ F?(x,yi 
J JAf 


Zo) d* dy 


(1) 


Implicit in this result are the following assumptions: the num- 
ber of structural modes within the frequency band of interest 
is large (i.e., approaches infinity), and the power spectrum of 
the wall response is slowly varying with respect to frequency 
relative to the rapidly varying acoustic modal transfer func- 
tion. As a result of these assumptions, the modal dynamics of 
the structure are effectively removed from the problem. In 
other words, the structural system is described in the AMA 
limit, whereas the acoustic system is taken in the CMA limit. 
Further, it is assumed that the modal damping ratio is small 
(<$1), which removes the coupling between acoustic modes 
(see Ref. 9 for justification). 

To obtain the AMA result for the acoustic cavity, further 
assume that the acoustic generalized mass squared (M? ) 2 , the 
frequency of the acoustic mode cubed (w?) 3 , and the acoustic 
damping (tf) do not vary rapidly with respect to modal num- 
ber r and can therefore be replaced by their values at the center 
frequency, (M* ) 2 , (o>?)\ and (tf ). Moreover, the expression 

Y i F}{x,y % z) F/(x,y,z 0 ) dx d y 
r J JA f 

is approximately equal to the average of F}{x,y,z) times 


e[U 2 

r JJ 


(x,y,z 0 ) dx d y 


as r — oo 


(i.e., a large number of acoustic modes). The previous expres- 
sion can be further simplified by 

E jj F r 2 (x t y,z 0 )dx dy = ^ A f (X^) A/ (Y^) A/ (Z^(z 0 )) A/ 

where the right-hand side reduces to 

(F?) 


A/AN a 


(Zt) 


where <Z C 2 ) = (F?)/(F C 2 ) A/ . Here, <F ? 2 ) is a volume average, 
and (F}) A/ is an average over the vibrating structural wall 
area. Then, 


(PnCl) 




**»{u t )A}{F*)ZF*lx,y,z) 


( 2 ) 


This result is the AMA representation for the acoustic cavity, 
which assumes that there are both a large number of acoustic 
modes and a large number of structural modes. It is valid at 
any point j r,y,z inside the acoustic cavity and is not a spatial 
average. The spatial average of Eq. (2) is 

<ft 2 ) ~ l( dfY m 

< p 0 cl ) 2 = 4 V V) Am* («? ) 2 <Z 2 > 

This expression is referred to as the spatially averaged AMA 
result ((AMA)) and is used as the denominator throughout the 
analysis. 

The effect of the source position is included in the term 
(Z 2 ), which is a ratio of the spatially averaged modal func- 
tions over the volume to that over the flexible wall area. Except 
for this term, the expression in Eq. (3) is independent of the 
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cavity geometry. The <Z C 2 ) ratio will always be of order 1. 
Specifically, for a rectangular geometry when the flexible por- 
tion is well removed from a corner or an edge, or for a cylin- 
drical cavity with a vibrating portion on one of the ends (but 
away from the edge), (Z 2 ) = Vi. The usual AM A result as- 
sumes this value of (Z c 2 ); however, the cases where the flexible 
portion is placed at a corner or an edge are easily accounted for 
in the (Z c 2 ) term and result in an additional multiplicative 
factor of Va or l A, respectively. 

In order to separate the effects of the response position 
in*>de the cavity from the position of the flexible portion Aj 
of the wall (source position), two ratios were used, (CM A}/ 
(AMA) and CMA Joca |/(AMA), Several combinations of CMA 
and AMA, spatially averaged and local results, were possible. 
The motivation for the choice of these particular ratios is 
presented in the Discussion section of this paper. The first 
ratio, comparing the spatial averages computed by the two 
methods, is used to demonstrate the convergence of the CMA 
results toward AMA as a function of frequency. The spatially 
averaged results are also used to demonstrate the overall effect 
of excitation position on sound pressure level. The second 
ratio is used to study local behavior at particular points inside 
v * acoustic cavity. This ratio is therefore used to demonstrate 
the effect of response position independent of source position, 
as well as the combined effects of source position and local 
response position. 

The ratio of the spatial average of CMA to the spatial aver- 
age of AMA is 



Flexible 

Vibrating 

Portion 


Fig. 2 Rectangular acoustic cavity with a portion of one wall flexible 
and vibrating. 


Equations (1-5) hold for any cavity geometry. 

Rectangular Cavity 

Dowell et al. 10 have shown that the acoustic pressure for a 
rectangular cavity with a flexible wall (all others rigid) can be 
described by the well-known rigid wall expansion or hard box 
modes for the structure: 


CMA _ y 
(AMA) " 



(x,y,z 0 )dx dy(u)?) 3 (Z?} 


A N a Aj 


(4) 


T'- .-, expression is derived in Ref. 4 and can be obtained from 
E (1) and (3). Equation (4) was used to assess the effect of 
excitation size and location due to area change and position of 
the vibrating portion of the wall. 

The separate effect of response position was studied by tak- 
ing the ratio of the local mean square pressure predicted with 
CMA to the spatially averaged mean square pressure predicted 
by AMA: 


(CMA) 

(AMA) 


(<^) 3 (Z C 2 ) 

AN A A f 


F;(x,y,z) F; l (x,y,Zo) dx dy 


(M r A ) 2 (u?y 


(5) 


Y 4 


FLEXIBLE 



Fig. 1 Flexible vibrating portion of one wall of the cavity. 


In this analysis, the flexible portion of the structural wall is 
allowed to vary both in size and position. Therefore, the inte- 
gral 



(x t y,z 0 )dx dy 


in Eqs. (4) and (5) becomes 



where r xt r y , and r z are modal indices, and x W q> y Wo , x w , y w , 
L xt and L y are defined in Fig. 1. This integral can then be 
evaluated analytically in terms of the parameters x Wq , y Wg and 
a w , b „ . In the previous work by Kubota et al., 4 which assumed 
that the entire wall was flexible, this integral over the flexible 
area A/ divided by Aj evaluated simply to Va , which, in gen- 
eral, is not true here since in most cases only a portion of the 
wall is flexible. 


Numerical Study 

Description 

For the numerical study, a rectangular cavity was considered 
(Fig. 2). One of its walls, or a portion thereof, was assumed to 
vibrate in an infinite number of structural modes. The wall was 
driven with banded white noise. The effects of varying both 
the size and position of the vibrating portion of the wall were 
studied. The size of the wall varied from full wall (100% wall 
area) down to a point (0.004% wall area). 

The quantities used in the study were ratios of mean square 
pressure predicted by CMA to that predicted by AMA. Ini- 
tially, spatial averages using both CMA and AMA were taken 
in order to avoid introducing the response location within the 
cavity as an additional parameter. Later, the local acoustic 
pressure responses at corner points, edge points, points on the 
face, and points in the interior were considered for these cases. 
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center frequency 


Fig. 3 (CMA)/(AMA) vs center frequency, 200-Hz fixed band- 
width, flexible area shrinks to center. 



center frequency 


Fig. 4 (CMA) /(AMA) vs center frequency, 400-Hz fixed band- 
width, flexible area shrinks to center. 

Results 

Spatially A veraged Case 

The ratio of the spatial average of CM A to the spatial aver- 
age of AMA [Eq. (4)] for the case where the oscillating portion 
of the wall converges to a center point is shown in Figs. 3-5. 
Each figure shows the spatially averaged CMA to spatially 
averaged AMA ratio for a different frequency bandwidth 
(200, 400, and 600 Hz) as a function of center frequency. The 
bandwidth Ao> is defined as Aoj — — oj m j n , and the center 

frequency oj c as u* = (a> ma x*a> m j n ), where ow and aw are the 
maximum and minimum frequencies of the frequency interval. 
All acoustic modes are assumed to have the same modal criti- 
cal damping ratio f 1). 

As can be seen from Figs. 3-5, all results approach unity as 
the center frequency becomes large. The larger bandwidths 
yield smoother curves, and the smaller bandwidths approach 
the asymptote slightly more rapidly. Kubota et al. 4 found sim- 
ilar trends regarding bandwidth effects in the related work. 
However, what was not expected was that departure from the 
entire wall oscillating caused litle change in the spatially aver- 
aged CMA/AMA ratio for the cavity. This was most likely due 
to the fact that the oscillating plate was centered about the 
midpoint on the wall and that all modes are symmetric or 
antisymmetric about that point. 


In Fig. 6, the result of the spatially averaged CMA/AMA 
ratio [Eq. (4)] for the oscillating plate of variable area and 
convergence to a point in the corner are shown. Again, plots 
that correspond to different frequency bandwidths yield sim- 
ilar results, and so only the 400-Hz fixed bandwidth plot is 
shown. In this case, there is a family of curves that approaches 
unity as center frequency (and therefore, number of modes) 
increases, as expected. However, the asymptote is approached 
from above rather than below for all plates smaller than the 
quarter wall. The quarter-wall case is equivalent, in terms of 
the CMA/AMA ratio, to the full wall due to symmetry. The 
cases where the plate is larger than a quarter panel approach 
from below as did the center point cases. For those cases in 
which the vibrating wall is smaller than a quarter panel, not 
only does the curve approach the asymptote from above, but 
as the oscillating portion of the wall better approximates a 
point, the peak of the curve approaches 4 and is slower to drop 
off to the asymptotic limit of 1. This region of elevated sound 
pressure level is similar to the intensification zones for struc- 
tural systems discussed in Crandall, 6 Itao and Crandall, 7 and 
Crandall and Kulvets 8 and for structural-acoustic systems in 
Kubota et al. 4 However, the intensification is due to excitation 
location rather than response location. 

It may first apear as if the AMA result does not account for 
this intensification. However, it was assumed that the excita- 



center frequency 


Fig. 5 (CMA) /(AMA) vs center frequency, 600-Hz fixed band- 
width, flexible area shrinks to center. 



center frequency 


Fig. 6 (CMA/(AMA) vs center frequency, 400-Hz fixed bandwidth, 
flexible area shrinks to corner. 
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Four Particular Points of Interest * denoted by A, B, C, D 



Flexible 

Vibrating 

Portion 


Fig. 7 Points at which the local response was predicted. 


CORNER PT (0, 0. 0) AS PLATE CONVERGES TO CENTER 


CMA/AMA 



8a Local CM A to (AM A) ratio vs center frequency for point A t 
flexible area shrinks to center. 


CORNER PT (0. 0. 0) AS PLATE CONVERGES TO CORNER 


CMA/AMA 



Fig. 8b Local CMA to (AM A) ratio vs center frequency for point A , 
flexible area shrinks to comer. 

tion occurs away from a corner or an edge of the wall in 
computing the integral over the flexible area in the AMA ex- 
pression. The effect of excitation intensification can be in- 
cluded in the AMA expression and, for corner excitation, 
would provide a factor of 4 increase. In addition to studying 
the effects of varying size and position of the oscillating por- 
tion of the wall in a spatially averaged sense, the local response 
was also calculated. 


Local Response 

The effect of position in the cavity was determined by con- 
sidering the local response. This was again done for both con- 


vergence of the vibrating plate from a full vibrating wall down 
to a vibrating point in the center and convergence to a point in 
a corner of the wall. The results are presented as a ratio of the 
mean square pressure evaluated at particular points in the 
cavity using CMA to the spatially averaged mean square pres- 
sure predicted by AMA. This ratio was computed for many 
center frequencies at a constant bandwidth of 400 Hz. There 
was no need to vary the bandwidth since Figs. 3-5 show little 
variation with bandwidth. The result was computed for vari- 
ous center frequencies in 200-Hz bandwidth increments up to 
a center frequency of 6000 Hz. Beyond 6000 Hz, 1000-Hz 
bandwidth increments were taken up to 11,000 Hz. This pro- 
duces a smoother looking curve beyond 6000 Hz due to the 
larger frequency bandwidth increments. 

Initially, four special response points were considered (see 
Fig. 7): the corner point (0,0,0), the midpoint of the flexible 
wall, the midpoint of the entire cavity, and a point on the wall 
along the centerline (1 .8, 1 .5 ,0). For these four points, the ratio 
of CMA to spatially averaged AMA was plotted as a function 
of center frequency in Figs. 8-11. 

Taking the corner point first (point A in Fig. 7), Fig. 8a 
shows that, as the plate converges to the center of the wall, the 
response of the corner approaches a pseudoasymptote of 8, 
whereas for convergence of the plate to the corner (Fig. 8b), 
this same point has a pseudoasymptote of 32, a factor of 4 
greater. The idea of a pseudoasymptote will be discussed in the 
following section. The factor of 4 comes from the excitation 

FLEX HALL MID PT AS PLATE CONVERGES TO CENTER 


CMA/AMA 



Fig. 9a Local CMA to (AMA) ratio vs center frequency for point B, 
flexible area shrinks to center. 


FLEX HALL MID PT AS PLATE CONVERGES TO CORNER 


CMA/AMA 



Fig. 9b Local CMA to (AMA) ratio vs center frequency for point B, 
flexible area shrinks to corner. 
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PT (1 .8, 1.5. 0.) ON WALL — CENTER CONVERGENCE 


CMA/AMA 



CENTER FREQUENCY 


Fig. 10 Local CMA to (AMA) ratio vs center frequency for point C, 
flexible area shrinks to center. 


CAVITY MID POINT AS PLATE CONVERGES TO CENTER 


CMA/AMA 



Fig. 11 Local CMA to (AMA) ratio vs center frequency for point D, 
flexible area shrinks to center. 


location (i.e., the location of the point source rather than the 
point of evaluation being in the corner), as was seen in the 
spatially averaged cases. Recall that when a point source was 
in the comer the spatial average was 4 times greater than when 
the point source was in the center of the wall. 

Figures 9a and 9b show that, at the midpoint of the flexible 
wall (point B in Fig. 7), both types of convergence yield a 
pseudoasymptote of 8. Since this point is on the wall, its ex- 
pected pseudoasymptote is 2. However, when the excitation is 
in the corner (Fig. 9b), this is increased by a factor of 4, hence 
the value 8. Five curves are actually plotted; only the curves 
representing the smallest plate areas (0.01 and 0.0049b) deviate 
significantly from the others. On the other hand, when the 
excitation location and the response location are at the center 
of the wall (Fig. 9a), the response there is also increased by a 
factor of 4. This phenomenon is similar to the intensification 
observed by Kubota and Dowell 3 in experiments where point 
loads are applied to a rectangular plate. They found that, in 
the limit of a large number of responding modes on the plate, 
the response of the plate as measured by the plate acceleration 
at various places on the plate surface was nearly uniform. The 
accelerations were significantly higher near the application 
points of the point loads. Because of symmetry, the effect of 


a point sound source acting at the center of the wall is to divide 
the rectangular cavity into quadrants (defined by cutting per- 
pendicular planes through the point source). In the newly de- 
fined subcavities, this point where the source is located is now 
a comer point. The response at a corner point is eight times 
greater than the interior region, and so the pseudoasymptote 
of 8 is appropriate for this case. 

Kubota and Dowell 3 also found hot lines running perpen- 
dicularly through the point force. To test for these in this 
acoustic cavity analysis, a point along one of the anticipated 
hot plates was studied (point C in Fig. 7). Since this point is on 
the face, it is expected to have a pseudoasymptote of 2.0, 
which is indeed the case if the full wall is moving (dotted line 
in Fig. 10). However, in the case of center convergence (i.e., 
analogous to a point load acting at the center of the structural 
wall), this point lies on a hot plane, and Fig. 10 shows a 
pseudoasymptote of 4. Assuming that the hot planes divide the 
cavity into subcavities, this point is on an edge of a subcavity. 
Therefore, the value of 4 is appropriate since the response at 
an edge is four times greater than the interior. When the oscil- 
lating portion of the wall converges to the corner, the pseudo- 
asymptote is 8. This is consistent with previous findings for 


ALONG THE EDGE - CENTER CONVERGENCE - {8000 HZ) 

CMA/AMA 



Z-0ISTANCE IN FEET 


Fig. 12 Local CMA to (AMA) ratio vs distance along an edge of 
the cavity away from a comer, with a point sound source in the center 
of the wall. 


RADIALLY INTO CAVITY - CENTER CONVERGENCE 

CMA/AMA 



RADIAL DISTANCE FT. 


Fig. 13 Local CMA to (AMA) ratio vs distance radially into the 
cavity away from a corner, with a point sound source in the center of 
the wall. 
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C ornef convergence since it represents a factor of 4 increase 
1 f over the expected value for a point on a wall. 

@ Response at the midpoint of the entire cavity (point D in 
Fig. 7) was also considered. At an interior point such as this, 

< the expected asymptote is 1.0. However, when the plate con- 
verges to the center of the wall (Fig. 11), this point lies on a 
plane of symmetry resulting in a factor of 4 increase. Similar 
to the previous case, this point now lies on an edge of a newly 
defined subcavity. The edge point response is four times 
greater than the interior, and so the center convergence case is 
consistent. When the plate converges to a corner of the wall, 
% i^ain a fourfold increase is expected, due to the previous re- 
■Jt- suits for the spatially averaged cases, and the result is a pseu- 
• r doasymptote of 4 (see curve E in Fig. 11). 

Points A-D were studied as the plate size was allowed to 
■ vary, for the two cases (sound source shrinking from a full 
vJ wall to a point sound source in a corner of the wall, and to a 

; > point in the center of the wall) as a function of center fre- 

quency (for a Fixed bandwidth). Next, the plate size was Fixed 
* at 0.004% of the wall area, corresponding to a vibrating point. 
The center frequency was Fixed at a value at which the pseu- 
doasymptotes had previously been reached (8000 Hz), and the 
»; and width was Fixed at 400 Hz. The distance into the cavity 
i . . ;m the vibrating wall was varied. Sample plots of local CM A 
to spatially averaged AMA vs distance into the cavity are 
shown in Figs. 12 and 13. 

In Fig. 12, the sound source (vibrating point) is located in 
the center of the wall, and the response is plotted along an 
edge. The peak response in the corner is 8. Moving away from 
the comer, the response then oscillates before approaching the 
asymptote for an edge, which is 4.0. The data are symmetric in 
the z direction, as can be shown analytically. Therefore, only 
half of the edge length is plotted (3.5 out of 7.0 ft). The 
purpose of this plot is to show that there is a region where the 
transition is made from the increased (intensification) value to 
ir s .\ nearly uniform value farther away. Further studies are 
n^zled to examine this transition zone and the parameters that 
determine its thickness. 

Figure 13 is a plot of the response in a radial direction away 
from the corner of the cavity for the sound source in the cen- 
ter of the wall. The radial direction is defined by the line 
i x -y — z , and the radial distance is equal to the square root of 
(x 2 +y 2 + z 2 ). Since the point source is in the center of the wall, 
hot lines exist that run down the center of the cavity. Because 
j of these hot lines, which redeFine new effective boundary 
points, the cavity interior is no longer uniform, as is shown in 
Fig. 13. After the radial distance of 1.0, the response begins to 
in ; ease, and approaches a value of 2.0, as if there were a wall 
Oi ;ice there. This is not a physical boundary created by the 
cavity geometry, but rather an artiFicial boundary created by 
the point source. If the same radial trajectory is taken with the 
point source located in the corner, the response of the interior 
is uniform. The peak value in the comer is 32; the response 
then oscillates and eventually approaches a uniform interior 
value of 4.0. This plot shows the theoretical existence of hot 
planes or surfaces inside the cavity that are created by a point 
source acting at the center of the wall. 

In summary, for a vibrating point at the center of the wall, 
the asymptotic limit for points that do not lie on hot planes is 
1.0 for interior points, 2.0 for points on a face, 4.0 for points 
on edge, and 8.0 for comer points. Also, the comer conver- 
gence cases yield the same relationships between locations but 
the magnitudes are increased by a factor of 4. Hot planes can 
be thought of as dividing the cavity into subcavities or quad- 
rants along lines of symmetry. Each subcavity then produces 
its own corner, edge, and face points, with the respective 
sound pressure levels relative to each other, thereby redefining 
effective boundary points. The intensiFication results that were 
obtained here were through the use of CMA; however, AMA 
is capable of predicting intensification as well, as is pointed out 
in the next section. In subsequent work, 11 intensification zones 
are studied using an extension of the AMA method. 


Discussion 

The term in the CMA/AMA equation [Eq. (4) or (5)] that is 
affected by changing the flexible area size and location of the 
flexible portion is 

F}(x,y,z 0 )dx dy/Aj- 
A i 

This can be thought of as a spatial average of the acoustic 
modal function in two dimensions ( x and y ). The expected 
result would be V * , unless the argument of one or both cosine 
functions (in F r ) is always zero or a multiple of tt. When the 
plate converges to a point in the corner, the x and y values are 
essentially zero; therefore, the value for the cosine is equal 
to 1 and the spatial average would then be 1.0 rather than !4. 
Mathematically, the effect of shrinking the area down to a 
point in the corner yields a fourfold increase in the relative 
mean square pressure locally or spatially averaged. 

However, one can imagine driving the frequency up so high 
that the approximated corner point no longer behaves like a 
point compared to an acoustic wavelength. This occurs when 
the acoustic wavelength becomes so small that the fixed phys- 
ical size of the sound source is much larger than the wavelength 
of sound at that frequency. It is for this reason that the term 
pseudoasymptote has been used rather than asymptote. As the 
center frequency becomes sufficiently large (i.e., wavelength 
sufficiently small), the true asymptote will always be 1.0 for 
the spatially averaged CMA/AMA ratio. 

Throughout this study, the numerical results have been 
presented either as the ratio of spatially averaged CMA to 
spatially averaged AMA, (CMA) /(AMA), or for the ratio of 
local CMA response to spatially averaged AMA response. 
Other possible ratios could have been used. An alternative 
choice could have been the ratio of local CMA to local AMA. 
However, this ratio would not show any variation between 
points on a boundary and points in the interior. It would 
approach 1.0 at high center frequencies for all points if the 
excitation were in the center, or 4.0 if the excitation were in the 
corner, as is shown in the following analysis. 

The local AMA result is [Eq. (2)] 

p 2 A f AN A {F 2 ) ^ F?(x,y,z) 

ip 0 c 2 ) 2 = 4 (M A ) 2 (u a ) 3 fc (Z 2 ) r AN A 

whereas the spatially averaged AMA result (derived in Ref, 4) 
is 

< p 2 ) _ r A N A (aX 2 Aoj<Mcj c ) 

(p 0 c 2 ) 2 % AuA VJ (atf) 3 tf (Z 2 ) 

Therefore, the ratio of local AMA to spatially averaged AMA 
is 

(AMA) local /(AMA) jpltW average “ \jjF}(x,y,z)/&N A ~^ / (F}) 
The numerator 

[E/vW.zyAN-*] 

is equal to V% when x t y , and z are not zero or L Xi L yy L z . 
It is equal to V* when one of the values of x t y , or z is equal 
to 0 or the length of the cavity in the appropriate direction, 
which is true on any face. Similarly, the numerator becomes Vi 
for an edge point and 1 for a corner point. 

The spatially averaged acoustic modal function evaluated 
at the center frequency (F 2 ), which comprises the denomina- 
tor of the (AMA) l0Cll /(AMA) jpltial avenge ratio, is always equal 
to Vt. Therefore, the (AMA)iocai/ (AMA) S p gt i 4 ] aV erage is summa- 
rized in Table 1. 
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Table 1 Local to spatial 
average-AMA ratios 


Location 

AM A local/ (AMA) 

Corner 

(l)/(l/8) = 8 

Edge 

(l/2)/(l/8) = 4 

Face 

(l/4)/(l/8) sa 2 

Interior 

(l/8)/(I/8)= 1 


Recall that the local CMA to spatially averaged AMA ratio 
for the center convergence case yields pseudoasymptotes of 
8 in the corner, 4 on an edge, 2 on a face, and 1 in the interior. 
These ratios are four times greater if the excitation is in a 
corner. 

This indicates that, for a large number of modes, the AMA 
results agree locally with the exact (CMA) results predicted 
when the oscillating wall is a full wall or converging toward the 
center. For corner convergence, the difference is a factor of 4. 
In deriving the AMA result used in this study, it was assumed 
that the excitation occurs at a location other than in a comer 
or on an edge, which accounts for the factor of 4 difference 
between center excitation and corner excitation, as explained 
previously. It is possible to incorporate the excitation location 
effect into the AMA result, if desired. 

Conclusions 

An AMA approach has been developed and applied to a 
coupled structural-acoustic problem. It is broadly applicable 
to any linear dynamic system at high frequency regardless of 
geometry, assuming that the damping is small. Further work is 
needed to develop an AMA result that is applicable for cases 
involving large damping. 

It is an extremely flexible approach and can be developed in 
accord with the nature of the system under study through 
inclusion of a series of simplifying assumptions. This tech- 
nique can thereby bridge the gap between CMA and SEA in 
terms of computational requirements and predictive capabil- 
ity. In that AMA is developed from CMA, it retains the capa- 
bility to predict spatial variations (intensification) in sound 
pressure levels or other relevant responses, something of which 
SEA is not capable. On the other hand, simplifications arising 
from the nature of the forces and the number of structural and 
acoustic modes involved result in a process that does not re- 
quire individual modal characteristics. This greatly reduces the 
number of calculations required relative to CMA. 

A rectangular acoustic cavity, with Five rigid walls, was 
chosen to investigate the capabilities of AMA. Spatial averages 
and local behavior for sound pressure levels were calculated 
for a number of cases involving the location and size of the 
sound source on the wall. For the spatially averaged cases, 
intensification due to source location was observed. In partic- 
ular, when a point sound source was located in the corner as 
opposed to the center of a wall, the spatially averaged sound 
pressure ratio was increased by a factor of 4. 


In addition to the spatial average, the local response was 
also calculated. The local response of the cavity interior is 
nearly uniform, with the exception of points on the structural 
boundary (walls, edges, and corners), when one entire wall of 
the rectangular cavity is vibrating. However, when only a por- 
tion of one wall vibrates, and particularly when this portion 
approaches a vibrating point, there are further exceptions. 
Perpendicular planes (hot planes) that run through the vibrat- 
ing point were found to divide the cavity into subcavities, 
which have new corners, edges, and walls. New subcavity cor- 
ners, edges, and walls exhibit the same relative increase as the 
original comers, edges and walls, which is eight, four, and two 
times greater than the interior, respectively. 
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Nomenclature 

area 

speed of sound 

cavity acoustic modal function 
: frequency 
= wave number 
= cavity dimension 
= generalized mass 
= modal index 
= pressure 

= real part of impedance 
= volume 
= displacement 

= spatial position coordinates 
= imaginary part of impedance 
= impedance at the boundary 
= number of acoustic modes 
= damping ratio 
= density 

= power spectrum 
= frequency 

= spatially averaged quantity 


Subscripts 

b = bandwidth 

c = center frequency 

/ = flexible 

o = reference value 

r = acoustic modal index 


Superscripts 

= acoustic 


A 


= time derivative 
= (overbar on pressure) rms 
= (overbar on impedance quantities) 
nondimensionalized by pc 
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Introduction 

I N structural-acoustic coupled systems, such as the interior 
of an automobile or an aircraft fuselage, an acoustic fie 
is created by the structural vibration of a wall or walls of th 
enclosure. This problem has been studied previously using 
asymptotic modal analysis (AMA) on a rectangular acoustic 
cavity u The cavity was assumed to be entirely rigid except for 
a vibrating portion of one of its walls The vibrating ; portion 
contained a large number of structural modes, which m turn 
generated a large number of acoustic modes in the interior 
space. Acoustic theory predicts, and the previous numerical 
work has shown, that there are intensification zones in the 
acoustic field near the cavity boundary and an otherwise uni- 
form response in the interior region. This is due to the lack of 
spatial correlation in the interior and the imposition of spatia: 
correlation at the boundaries. For a rectangular acoustic ca - 
ity, it is well known that the mean square pressures are lg , 
four, and two times the uniform interior pressure levels at the 
corners, edges, and faces, respectively. 3 ' 

In designing acoustic spaces, allowances must be made for 
these intensification zones. Therefore, it is important to deter- 
mine the characteristic distance over which the response levels 
change from their peak values at the boundary to the uniform 
interior level. If the distance from a boundary is nondimen- 
sionalized by the center frequency wave number k c , then the 
intensification zone can be described by the nond.mensional- 
ized spatial variable k<x and the bandwidth to center fr 
quency ratio/*//,. Significantly, the parametric dependence is 
independent of cavity dimensions. The spatial variation of 
mean-square pressure is also independent of cavity dimen- 
sions. The spatial variation of mean-square pressure is a so 
independent of the size of vibrating or absorptive surfaces on 
the cavity wall, provided that the surfaces are very large com- 
pared to an acoustic wavelength (high-frequency limit) and 
provided that the damping is small. 

The dependence on bandwidth can be separated mathemati- 
cally from the dependence on center frequency through a 
Taylor series expansion in f b /f c taken about ^e center re- 
quency wave number. It is found that the bandwi t e P 1 

dence is a higher order effect. The resulting mathemauca 
expression for mean-square pressure consists of a term that s 
dependent only on the center frequency, plus 'erms that co 
tain the bandwidth to center frequency ratio; but these term 
are of order 0[(/c//c) 2 ] ar >b higher. 


1 
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Using AMA, the spatial variation of mean-square pressure 
was found to be independent of cavity dimensions. This sug- 
gests that the problem can also be analyzed as a local problem 
consisting of an infinite number of oblique acoustic waves 
incident upon a rigid surface from all possible angles. The 
solution to that problem is equivalent to that derived from 
AMA. 

In the same manner, an expression can be obtained for the 
spatial variation of mean-square pressure, which describes the 
intensification zone near a boundary for a case where the 
incident wall is nonrigid (i.e., finite impedance cases). This 
result can be written in terms of a hard wall component plus 
an absorption correction term, which contains the impedance 
information. 

The purpose of this research is to determine the structure of 
these intensification zones and their relation to the entire 
acoustic field. It is hoped that this work will lead to better 
design and analysis of acoustic spaces by providing a more 
thorough understanding of intensification zones for both rigid 
and absorptive walls. 


Background 

Intensification zones for structures have been studied exten- 
sively by Crandall 6 and others for random vibrations of plates 
(see, e.g., Refs. 6-8). They have found that the response of the 
plate at high frequency is relatively uniform, with exceptions 
occurring at the point of application of a force and/or at the 
boundaries, depending on loading conditions and boundary 
conditions. Kubota and Dowell 9 found similar intensification 
phenomena for a rectangular plate excited by random vibra- 
tion of a point force. They used AMA to predict the response 
of the plate, and they verified their results with experimental 
data. 

The study of intensification zones in acoustic spaces was 
first begun by Waterhouse 10 in the 1950s. His pioneering work 
on interference patterns in reverberant sound fields was fol- 
lowed by further research by Chu, 11,12 Tohyama and Suzuki, 13 
and others 14,15 who were specifically interested in the place- 
ment of microphones in reverberation chambers. The results 
presented here for rigid wall boundaries agree with the earlier 
work, although the derivations are different. References 10-15 
were restricted to rigid walls only. Later work by Waterhouse 
and Cook 16 included the treatment of pressure release (p - 0) 
boundary conditions. 


Parameterization of Intensification Zone 
Asymptotic Modal Analysis Approach 

A rectangular acoustic cavity with five rigid walls and one 
vibrating wall has been studied previously using AMA. A 
diagram of this problem is shown in Fig. 1. The vibrating wall 
was driven by white noise such that a large number of struc- 
tural modes was excited in a particular bandwidth. It was 
assumed that a large number (approaching infinity) of acous- 
tic modes was present in the interior acoustic cavity, and that 
these modes were temporally uncorrelated. Although it is not 
assumed in the AMA derivation, a consequence of the previ- 
ous assumptions is that the acoustic modes are spatially uncor- 
related as well. From Refs. 1 and 2, the expression for the 
mean-square pressure for a finite number of acoustic modes 
and an infinite number of structural modes is 


P 2 

(fioCo) 2 




F?(x.y>z) 


x 



( *>y,z 0 ) dx d y 


(1) 


As in AMA, M? and can be treated as constant in a band 
and can therefore be evaluated at the center frequency. Thus, 
for example, becomes A/^, etc. Normalizing p 2 by the 
spatially averaged p 2 yields the following expression for the 


nondimensionalized mean-square pressure in the rigid cavity 
with one entire wall flexible and vibrating: 

pi 

T 

where the quantities in () are spatially averaged values. This is 
the expression that describes the intensification zone near a 
boundary; it is a ratio of the local mean-square pressure to the 
spatially averaged (or uniform interior value) mean-square 
pressure. At interior points, this ratio will asymptotically (for 
r — oo) approach 1.0, whereas at the boundaries, the ratio will 
be 2.0, 4.0, or 8.0 for walls, edges, or corners, respectively. 
There are several possible transitions to study: from a corner 
to an edge, a wall, or the interior; from an edge to a wall or the 
interior; and from a wall to the interior. For each case, there 
are one-, two-, and three-dimensional modes to consider. 

For the rectangular acoustic cavity, Dowell et al. 17 have 
shown that hard-box modes can be used to describe the inte- 
rior acoustic field, even when all six walls are not rigid. There- 
fore, 



One-Dimensional 

Considering one-dimensional modes only, F r = cos ( mrx / 
L x ) and cj? = mrc/L x . Assuming there are a large number of 
acoustic modes, which is consistent with the AMA approach, 
the discrete variable n can be treated as a continuous variable. 
Equation (2) becomes, 




Fig. 1 Rectangular acoustic cavity that was studied. 
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1-D Intensification 
Summation Vs. Integration 



Fie. 2 Nondimensionalized mean-square pressure vs k c x computed 
three ways: as an integration in wave-number space, as a summation 
over 15 modes, and as a summation over 567 modes. 


where f c is the center frequency, f b the frequency bandwidth, 
and f ( and/ u are defined as the lower and upper frequencies, 
respectivel y, of the frequency interval as follows: f u -ft = fb> 
and f c = JjTlt. The arguments of the cosines and the cosine 
integral function Ci are defined by 0 f - 2*k c x*ft/f c , 0u- 
2 *k c x*f u /f C i where the nondimensional distance k c x is equal 

'he ratios ///c and/ M // c can be obtained from f b /f c using 
these definitions. In fact, the entire expression (8) can be 
written in terms of f b /f e and k c x, where Ayr is the product of 
the wave number associated with the center frequency and the 
distance from the endpoint. Knowing the ratio of frequency 
bandwidth to center frequency (fb/fc) allows the pressure 
function in the transition zone to be plotted as a function of 
nondimensionalized distance M away from the endpoint. 
Note that the dimensions of the cavity do not appear in this 
result. 

If we assume that the frequency is approximately constant 
over the interval, and allow « to equal « c , the result would be 


IL 

{ P 2 ) 


1 1 

1 + 2 (kcX)(f b /f c ) 


{sin (0 U ) - sin (0f)) 


(5) 


This is a decaying function that beats with some predictable 
frequency, as is shown in Fig. 2. In Fig. 2, both the previous 
approximation (5) and the discrete sum are plotted. The sum- 
mation is performed for two different center frequencies, one 
in which 15 modes are present and the other in which 567 
modes are present. The summation expression requires an 
additional parameter that is dependent on cavity dimension. 
Although the AMA (integration) approximation is theoreti- 
cally valid when there are a large number of modes, here the 
ap;* oximation works well for relatively few modes. However, 
the agreement is not as good in the two- and three-dimensional 
cases. 


Two-Dimensional Case 

The transition zone for the two-dimensional case is slightly 
more complicated. In this case, F r » cos (nxx/L x ) cos (mry/ 


Ly) and (<^) J = c 2 [(mcx/L x ) 2 + (mry/L y ) 2 }. Assuming there 
are a sufficiently large number of modes, m and n can be 
treated as continuous variables and the summation can be 
replaced by integration. Converting to wave-number space 
( nrx/L x =k x , etc.) and transferring k x and k, into polar 
coordinates by the relations k x = k cos a and k y = k sin a, Eq. 
(2) can be written for the two-dimensional case as, 


Pi 

(P 2 ) 



[cos 2 (kx cos a)] [cos 2 (ky sin a)]/k 2 da d k 



da d k 


( 6 ) 


To show that this integral is dependent only on k<jc and fb/fc > 
let { = kx. Then ky = £ y/x and d£ = x dk. The nondimen- 
sionalized pressure ratio becomes 





(k u /kc)k c x 


[cos 2 (£ cos a)] [cos 2 (£y/x sin a )]/£2 da 




j* f j j 

— - d£ da 

<*/*r)*c* 4 £ 


1 


(7) 


In Eq. (7), the only independent parameters are k<x and the 
ratios k u /k c and k,/k c , which can be derived homf b // c - The 
ratio y/x is known from the desired location within the cavity. 
The geometry of the cavity itself (i.e., the cavity dimensions) 
is not an important parameter in determining the nondimen- 
sionalized pressure ratio of the cavity interior. 

If we assume that it is a constant (i.e., Ak /k < 1), this 
expression can be integrated explicity with the result that 

^sl+ J 0 (2k c x) + J 0 (2k c y) + J 0 (2k c x) ■ J 0 (2k c y) 

ip 1 ) 


+ ~ t t (- l)" Jm(2k c x)J lm i2k c y) 

X m ■ 1 n * 1 


^si n (m — n}x/ 2 sin (m 4- n)x/2) m ^ n 

2(m - n) 2(m + n ) 

(7t/4) m = n 


Notice in the limit as k<jc and k c y asymptotically approach 
infinity, the Bessel functions decay and the nondimensional- 
ized pressure ratio asymptotically approaches I, whereas the 
k c and k<y approach zero (i.e., near the corner), the higher- 
order Bessel functions decay, but Jq(0) = 1 and so the nondi- 
mensionalized pressure ratio asymptotically approaches 4. 

Three-Dimensional Case 

For the three-dimensional case, F r = cos (nxx/L x ) cos 
(rmry/L y ) cos (lxz/L : ) and u 2 = c 2 {(n*/L x ) 2 + (mx/L y ) 2 + 
(h/Lf) 1 ]. Substituting these expressions into Eq. (2), making 
the appropriate assumptions to allow the summation to be 
replaced by integration (large number of modes, etc.), and 
transforming from modal index space to wave-number space 
yields the following expression in spherical wave-number co- 
ordinates: 


8 

Pl = A 

(P 2 ) 


[cos 2 (^x sin 0 cos <t>) • cos \ky sin 0 sin <p) • cos 2 (A:z cos B)/k] sin 0 d0 d<f> dk 


J k,*,e 


(1 /k) sin 0 d0 d<f> dk 


k,*,e 


( 8 ) 
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1-D Intensification Curves 



Fig. 3 Nondimensionaiized mean-square pressure vs k c x one-dimen- 
sional case, three f b /f c ratios plotted: fb/fc = 0.005, 0.23, and 0,50. 


2-D Intensification Curves 


rate curves are beginning to converge (Fig. 4), whereas in the 
three-dimensional case (Fig. 5), convergence occurs more 
rapidly. This trend is due to the large number of two- and 
three-dimensional modes for a given center frequency at fixed 
bandwidth. 


Oblique-Wave Approach 

Alternatively, the problem of intensification near 
boundaries can be analyzed as a local problem. At the walls, 
for example, the problem can be modeled as an infinite num- 
ber of sound waves incident upon a rigid surface from in- 
finitely many directions. The expression for the pressure at a 
point in the x,y plane due to a single sound wave, as shown in 
Fig. 6, is p(x,y,t) = P where P = 2 A cos k x x t and 

A is the amplitude of the incident pressure wave. Assuming 
equal amplitude incident waves with random phase and ran- 
dom incidence angles (all equally probable), the summation of 
the mean-square pressures is 


£ P 2 = E 2 A 2 cos 2 (kx cos 6) 

A±j,8 A k,8 

where k x has been replaced by k cos 6. 

Assuming there are a large number of frequencies in a 
particular bandwidth, as in AMA, the mean-square pressure 
in the band is 



' k upper 
k lower 


cos 2 (kx cos 0) d k dd 



Fig. 4 Nondimenslonalized mean-square pressure vs k c x , two-dimen- 
sional case, three f b /f c ratios plotted: f b /f c = 0.005 , 0.23, and 0.50. 


Again, it can be shown that the nondimensional pressure 
ratio is independent of cavity geometry and is, in fact, only 
dependent on the parameters k ^ y/x, z/x, and f b /f c by per- 
forming a change of variables, as in the two-dimensional case. 
When the spatial variables are nondimensionalized by wave 
number, the transition zone can be plotted as a function of 
mean-square pressure vs position if the bandwidth to center 
frequency ratio is specified. 

As the number of dimensions increases, the ratio /<,// c plays 
a less significant role in the shape of the transition zone, as can 
be seen upon comparison of Figs. 3-5. In each of the three 
figures, curves corresponding to three different f b /f c ratios are 
overlayed. The three f b /f c ratios are f b /f c = 0.005 (which cor- 
responds to a narrow bandwidth at a high center frequency), 
f b /f c = 0.23 (corresponding to a X A octave band), and fj c - 
0.500 (which corresponds to an octave band). For the one-di- 
mensional case, Fig. 3 shows three distinct curves for the three 
different ratios. In the two-dimensional case, the three sepa- 


3-D Intensification Curves 



Fig. 5 Nondimenslonalized mean-square pressure vs k c x , three-di- 
mensional case, three f b /f c ratios plotted: f b /f c = 0.005, 0.23, and 
0.50. 




Fig. 6 Reflection from a rigid surface (left) or from an absorptive 
surface (right). 
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In order to compare this result to the AM A result, it must be 
nondimensionalized by its asymptotic far-field value (i.e., 
large c), which is 2* A 2 (k u -k i). The corresponding result 
derived from AMA is Eq. (6). However, Eq. (6) contains 
variation in both the y and * directions. The nondimensional- 
ized mean-square pressure ratio is equivalent to Eq. (6) for the 
two-dimensional AMA intensification from a wall providing 
that the y variable is held constant and the wave number is 
assumed large, such that the product ky approaches infinity. 
\ote that the 1 /k 2 terms in Eq. (6) are contained in the 
pressure amplitude P for the oblique waves. These results also 
agree with the earlier result of Waterhouse 10 which was 
derived using methods similar to the oblique-wave formula- 

U °Similar relations can be derived for the two- and three-di- 
mensional incident sound fields with reflections from one, 
two, or three rigid surfaces. The methodology is similar to that 
presented here, and the results are in agreement with the AMA 
results. 

Absorptive Wall Intensification 

Using the oblique-wave approach outlined earlier, the effect 
. i an absorptive wall can also be analyzed. Assume that an 
infinite number of oblique waves are incident upon an absorp- 
tive surface (of finite impedance, z b ) from infinitely many 
directions. Figure 6 shows one such pressure wave. The ex- 
pression for the pressure from one wave is 

f e -ikxx Y 

p = Pe^'-Kyy) = 2/1 COS kjc — = — -\e‘^'- kry) 

y l z b cos e + u 

and 

Z b = z b /p„c 

Rewriting, p = IP where P = \P\e i *". Then 
p „, |= IP I cos (a >t-k y y + <t>p). The mean-square pressure is 
equal to 

No. of wave* I P 1 2 

/> 2 = z -^r 

J - 1 2 

As in the previous analysis, assume that there are a large 
number of modes and replace the summation with integration. 
For equal amplitude waves, the expression for the mean- 
square pressure is 


Absorption Results: 2000 Hz Center Freq. 



Fig. 7 Comparison of theoretically computed and experimentally 
measured nondimensionalized mean-square pressure in dB vs k c x : 
2000-Hz center frequency, Vs octave band. 


easy to implement. This involves changing k c from a constant, 
evaluated at the center frequency, to a A:, a variable, and 
integrating the expression over the band Ak. 

Comparisons between theoretically calculated results and 
experimentally measured data are shown in Figs. 7-9 for three 
different center frequencies: 2000, 3250, and 4000 Hz, respec- 
tively. The impedance of the foam is different at each center 
frequency. At 2000 Hz, the normalized resistence (r b ) is 1.17 
and the normalized reactance is - 0.84. These imped- 
ance characteristics correspond to an equivalent random inci- 
dence absorption coefficient a of 0.85. At 3250 Hz, r b = 1 .44, 
*7 = 0.15, and the corresponding a = 0.95. At 4000 Hz, 
r b = 1.47, T b = 0.64, and the equivalent value for a is 0.91. 

The theoretical results are for the three-dimensional case 
with bandwidth included. These three examples are plotted for 
a l A octave bandwidth. The experimental data were taken as 
an extension of an experiment that was used to validate the 
asymptotic modal analysis method for a rectangular acoustic 


U 


p 2 = 2 \A H ) cos 2 (kcX cos 0) - 2 cos (J k c x cos 0) 


(F h cos 0 + 1) cos (kcX cos 0) - x b cos 0 sin ( k c x cos 0) 
(T b cos 0 + l) 2 + (x~ b cos 0) 2 


+ — 


( r b cos 0 + 1) ; 


i — ] 

) 2 + (x b cos 6) 2 ) 


dd 


Consistent with the assumptions of AMA, the wave number 
has been replaced by its value at the center frequency, assum- 
ing a narrow band with many modes. The first term in the 
ex:v ission is the same as in the rigid wall case. The far-field 
value (i.e., the limit for large k^) is 


P 2 = 2U 


1 2 ( r h cos 0 ' 

~ (F b cos 0 + l) 2 + (Xi cos 6)\ 


d0 


After dividing by the expression for the far-field mean-square 
pressure, the end result is an expression that is equivalent to 
the rigid wall mean-square pressure result, but with a set of 
terms added to include the effects of absorption. 

These expressions can be extended to three dimensions by 
integrating in spherical coordinates around both angular di- 
rections. If desired, inclusion of the bandwidth effect is also 


cavity with rigid walls. The apparatus and details of the exper- 
iment are described in Ref. 18. However, the original appara- 
tus was modified by placing an open-cell acoustic foam on one 
wall of the cavity. The experimental data were taken by mov- 
ing a microphone into the cavity from the foam-covered wall 
and recording the measured sound pressure levels. 

The theoretical curves show clear asymptotic behavior, 
whereas the experimental curves settle down less rapidly. The 
physical system in the experiment was a three-dimensional 
enclosure with rigid walls. Provided that the rigid walls were in 
the far field of the microphone, so that their intensification 
zones were not penetrated by the microphone, the experiment 
should be analogous to the theoretical model. Although this 
condition was satisfied, it appears from the experimental data 
that the interior of the cavity was not uniform. Previous work, 
both analytical 1 ’ 2 and experimental 18 have shown that, when 
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Absorption Results: 3250 Hz Center Freq. 



Fig. 8 Comparison of theoretically computed and experimentally 
measured nondimensionalized mean-square pressure in dB vs k c xi 
3250-Hz center frequency, Vi octave band. 


Absorption Results: 4000 Hz Center Freq. 



Fig. 9 Comparison of theoretically computed and experimentally 
measured nondimensionalized mean-square pressure in dB vs * c x: 
4000-Hz center frequency, Vs octave band. 


the damping in the cavity is small, the cavity interior is uni- 
form except in the intensification regions. Covering one wall 
with the absorptive foam apparently added enough damping 
to prevent the interior sound pressure levels from becoming 
uniform. 

The wall impedance dramatically changes the character of 
the intensification region. In fact, as Figs. 7-9 show, the 
sound pressure levels at the wall can be reduced well below the 
interior levels if the wall is sufficiently absorptive. Despite the 
discrepancies away from the wall, the curves show reasonable 
qualitative agreement between the analytical and experimental 
results near the wall, as is consistent with the local nature of 
the intensification process. 

Bandwidth Effect 

The key parameter in all of the previous studies was the 
bandwidth to center frequency ratio f b /f c - In order to isolate 


the contribution of the bandwidth from that of the center 
frequency, a Taylor series expansion was performed. For ex- 
ample, Eq. (6) for the two-dimensional waves in the intensifi- 
cation zone is a double integral in k and a. It can be written as 
a Taylor series expansion about the center frequency wave 
number k c as follows: 

rf 

F(k t a) d* da = F(k c , a) d k da 

k t JO J k ( 




dF(k,ct) 

dk 


(k — k c ) dk da 
kc 


[! f*“ 3 2 F(Ar,a) 

+ Jo J* r 


2! 


d/: da + - - • 


( 9 ) 


Using the arithmetic definition of center frequency, i.e., 
k c = ( k u + *j)/2, the integrals of odd powers of ( k-k c ) will 
equal zero, leaving only the integrals of even powers of band- 
width. Nondimensionalizing introduces a factor of bandwidth 
in the denominator. Therefore, the resulting function will 
consist of a term that is dependent only on center frequency 
plus higher-order terms that contain bandwidth squared and 
the higher even powers of bandwidth. 

An illustrative example, which can be solved in closed form, 
is the intensification in a corner approaching along an edge, 
where the motion of the wall is such that the displacement is 
relatively constant (white noise) in the bandwidth. For this 
case, the nondimensional mean-square pressure expression is 


(p 2 ) 



cos\kz cos 8) 


k sin 8 dd d <t> dAr 



k sin 6 dd d 4> dfc 


( 10 ) 


This equation differs from Eq. (8) by a factor of k 2 in the 
integrals, since Eq. (8) is for the case where the white noise 
assumption is applied to the acceleration of the moving wall, 
rather than its displacement. The z axis is taken as the edge 
along which the intensification is studied. This expression has 
the closed-form solution: 


p 2 _ 4 + 2 sin (2 k c z) sin [(k b /k c )k c z)] 
W)“ (k c z) 2 (k b /k c ) 


All center frequency subscripts here refer to the arithmetic 
center frequency. Alternatively, the Taylor series expansion 
procedure outlined in Eq. (9) and performed on Eq. 10 yields 


Pi 

ip 1 ) 


= 4 + 4 


sin(2£ c z) 
2 k c z 



(k c z) sin (2 k c z) 


+ higher order terms 


( 12 ) 


Expanding the sine ( k b /k c k c z) term in Eq. (11) in a Taylor 
series about the center frequency wave number produces an 
identical result. This approximate solution is only valid in the 
region where the sine function is closely approximated by a 
two-term Taylor series. For greater accuracy, the solution 
would have to include higher-order bandwidth to center fre- 
quency ratio terms. The two-term Taylor series for sine is only 
valid up to arguments approximately equal to 1 .25. Therefore, 
in this example, the combination ( k b /k c )(k c z ) must be ^ 1.25, 
or conservatively ^ 1.0, which means that k c z ^ 1 /(k b /k c ) 
for the approximate solution to be accurate. Typically, inten- 
sification occurs close to the wall, edge, or corner where the 





Fig. 11 Correction curve that is to be multiplied by f b /f c ratio 
squared: this curve by itself is also independent of bandwidth. 


Fig. 13 Comparison of exact solution vs the approximate solution 
consisting of CF curve + (/*// e ) 2 correction curve, for f b /f c = 0.500. 


values of k c x t k c y , or k c z are small. In fact, usually, k c x, k c y , 
and/or k c z £ 2x/3. Therefore, in most cases, the two-term 
solution should be sufficient. For a VS octave band, for exam- 
ple k b /k c is 0.23, which satisfies the requirement that k c z < 1/ 
0 k b /k c ). 

Since the bandwidth correction is second order in band- 
width, the mean-square pressure can be approximated by a 
si ..pie curve that is dependent on center frequency (CF) only 
for small bandwidths. This curve (the CF curve) for the cur- 
rent example problem is plotted in Fig. 10. It is a function of 
k c x alone, and is independent of bandwidth. The effects of 
increasing bandwidth can be taken into account by adding 
another function of center frequency that has been multiplied 
by the ratio (f b /f c ) 2 - This function, called the correction curve, 
is shown in Fig. 11 [before multiplication by (f b /fc) 2 ] for this 
example. 

In Fig. 12, the correction curve has been multiplied by 
(0.23) 2 and then added to the CF curve. The results are ex- 
pected to be accurate up to k c z around 4.0, as explained 


earlier. However, the exact solution and the approximate solu- 
tion agree with each other well beyond 4.0. In Fig. 13, a 
similar comparison between the approximate solution and the 
exact solution is made for f b /f c = 0.5. Here, good agreement is 
expected up to 2.0. Once again, the agreement actually extends 
beyond that predicted. From the plot, it is also easy to see that 
the approximate solution fails for large k c z. 

Conclusions 

In previous work, asymptotic modal analysis was used to 
predict interior sound pressure levels in a rectangular acoustic 
cavity. Local response peaks in sound pressure level were 
found at the boundaries; otherwise, the interior levels were 
nearly constant. In the present work, the spatial variation of 
sound pressure level between the boundary and the uniform 
interior was studied, i.e., the intensification zone. 

It was found that using such AMA techniques as treating 
the discrete summation as an integration, and evaluating cer- 
tain parameters at the center frequency, etc., provides accu- 
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rate results even when the number of modes is relatively few. 
In addition, the AMA formulation allows estimation of the 
mean-square pressure ratio as a function of nondimensional- 
ized distance into the cavity, without knowing the cavity 
dimensions. The family parameter for the curves of nondi- 
mensionalized mean-square pressure vs nondimensionalized 
distance (distance multiplied by the center frequency wave 
number) is the bandwidth to center frequency ratio /*// c . 
Plots of one-, two-, and three-dimensional intensification 
zones have shown that the dependence on fj c is less important 
as the number of modal dimensions increases, i.e., the effect 
is least important in the three-dimensional case. 

The bandwidth dependence was separated from the center 
frequency dependence through a Taylor series expansion 
about the center frequency wave number, which was per- 
formed on the expression for mean-square pressure. As a 
result, the mean-square pressure can be expressed as a func- 
tion of center frequency plus terms that are of order band- 
width to center frequency ratio squared, and higher (even 
powers). Therefore, if the bandwidth to center frequency ratio 
is sufficiently small, the expression can be simplified. If sev- 
eral bandwidth to center frequency ratios are to be considered, 
their effects can be added on as corrections that are of order 

(/i//c) 2 . 

The intensification zones were also analyzed as local prob- 
lems consisting of an infinite number of oblique incidence 
sound waves impinging upon a rigid or an absorptive surface. 
The result for the spatial average of the sound pressure levels 
in the rigid wall case was identical to that obtained from the 
AMA methods. The result for the absorptive wall was found 
to contain the rigid wall result as well as additional terms that 
contain the impedance information. The shape of the intensi- 
fication zone is strongly affected by the choice of wall 
impedance and can be remarkably different from that near a 
rigid wall. In particular, minimum rather than maximum lev- 
els may occur in this zone. The absorptive wall analytical 
result was compared with experimental results and found to 
agree qualitatively near the wall. 

The results of this work provide physical insight into the 
intensification zone behavior in acoustic cavities by providing 
analytical expressions that describe the sound field in terms of 
universal components that show the importance of various 
parameters. The analytical techniques that were used here 
provide efficient and accurate methods for predicting sound 
pressure levels in rooms with intensification zones. 
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Experimental Verification of the 
Asymptotic Modal Analysis 
Method as Applied to a 
Rectangular Acoustic Cavity 
Excited by Structural Vibration 

An experiment was performed on a rigid wall rectangular acoustic cavity driven by 
a flexible plate mounted in a quarter_ of one end wall and excited by white noise. 
The experiment was designed so that the assumptions of Asymptotic Modal Analysis 
{AM A ) were satisfied for certain bandwidths and center frequencies. Measurements 
of sound pressure levels at points along the boundaries and incrementally into the 
interior were taken. These were compared with the theoretical results predicted with 
AMA, and found to be in good agreement, particularly for moderate (1/3 octave) 
bandwidths and sufficiently high center frequencies. Sound pressure level measure- 
ments were also taken well in tot he cavity interior at various points along the 5 
totally rigid walls. The AMA theory, including boundary intensification effects, 
waw shown to be accurate provided the assumption of large number of acoustic 
modes is satisfied, and variables such as power spectra of the wall acceleration, 
frequency, and damping are slowly varying in the frequency bandwidth. 


Introduction 

Enclosed spaces such as the interiors of airplanes or auto- 
mobiles, in which the enclosure itself or portions of a wall 
vibrate, are most often subject to interior noise problems. This 
structural vibration induces an acoustic field in the interior 
space, and the wall flexibility allows external sound fields to 
be transmitted to the interior. Two methods which are tradi- 
tionally used to predict the interior noise levels are Classical 
Modal Analysis (CMA) and Statistical Energy Analysis (SEA). 
Another method, Asymptotic Modal Analysis (AMA) which 
contains aspects of both CMA and SEA, has recently been 
developed for structural-acoustic applications. 

Classical Modal Analysis involves consideration of each in- 
dividual mode of the system: the structural modes of the en- 
closure, the acoustic modes of the cavity, as well as the coupled 
structural-acoustic modes. Recent work by Pan and Bies (1990), 
studying the effect of fluid-structural coupling on sound waves 
in an enclosure is an example of a CMA-type approach. Other 
CMA-based methods include finite element analysis (for ex- 
ample, Sung and Nefske, 1984), and boundary element meth- 
ods (for example, Seybert and Cheng, 1987). Because individual 
modal contributions are considered, these methods are cum- 
bersome for systems with large numbers of modes. 
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The SEA method is often a good alternative method when 
there is a large number of modes involved, since only statistical 
quantities regarding the system need to be known. The SEA 
method gives accurate results, when the underlying stochastic 
assumptions are met in the frequency range considered. How- 
ever, because the results are given in terms of statistical quan- 
tities only (e.g., spatially-averaged mean-square pressures), 
detailed information cannot be obtained using SEA. The stand- 
ard reference on SEA is the textbook by Lyon (1975). 

A new method, Asymptotic Modal Analysis (AMA), which 
was developed by Dowell (1983), is a hybrid of the previous 
methods. It is systematically derived from Classical Modal 
Analysis, and therefore contains several levels of approxi- 
mation and generality. The accuracy and simplicity of the 
results obtained from AMA are adjusted by determining the 
level of approximation in the derivation. In this work, the 
accuracy and limitations of the AMA method for a structural- 
acoustic system are explored experimentally. 

Peretti and Dowell (1992a) analyzed the acoustic field inside 
a rectangular cavity where five walls were rigid and the sixth 
wall, or a portion of that wall, was vibrating. In their work, 
the Interior noise levels were studied through a comparison of 
the results obtained by Asymptotic Modal Analysis (AMA) 
and those obtained by Classical Modal Analysis (CMA). It 
was found that AMA, which in theory is valid in the limit 
when there is an infinite number ol modes responding, gave 


suprisingly good results when the number of modes was infinite 
and only moderately large. Perelti and Dowell (1992a) studied 
the case where the vibrating wall was excited with a white noise 
time history and assumed to respond with an infinite number 
of structural modes. Their theoretical study found that the 
sound pressure levels in the interior of the rectangular acoustic 
cavity were nearly constant, with the exceptions occurring at 
the boundaries. In particular, these boundary areas or "in- 
tensification zones” exhibited levels which were 8, 4, and 2 
times higher than the interior level for the corners, edges and 
walls, respectively. The intensification effect is caused by the 
existence of spatially correlated modes at the boundaries due 
to enforcement of the wall boundary conditions, while the 
modes in the interior are uncorrelated. , , 

Peretti and Dowell (1992b) also analyzed the transition be- 
tween the intensification at the boundaries and the uniform 
interior region. They found that the pressure field can be de- 
scribed in terms of a nondimensional parameter, which sca- 
lesthe distance into the cavity with the wavenumber 
corresponding to the center frequency, for different values of 
the bandwidth to center frequency ratio. This result is inde- 
pendent of the cavity dimensions for frequencies sufficiently 
high to achieve the AMA limit, and suggests that it is effectively 
a local solution near the boundary. In fact, similar intensifi- 
cation results were obtained for a uniform distribution of 
oblique waves incident on an infinite wall, which is the equiv- 
alent problem in the nearfield. 

The purpose of this work is to compare these previous an- 
alytical results with experimental data and to test the practical 
limitations of the theory. To this end, an experiment was per- 
formed in which a portion of one wall of an otherwise rigid 
rectangular cavity was excited structurally with a white noise 
time history. The experimented design was chosen to best match 
the assumptions which were made in the theoretical study of 
Peretti and Dowell (1992a). The principal findings of the earlier 
theoretical study which were further investigated by the present 
experimental study are: that the response in the interior of the 
cavity is uniform, that at the corners, edges and walls the sound 
pressure levels are 8, 4, and 2 times greater than the levels in 
the interior; and that the thickness of the intensification zone 
is independent of the cavity dimensions. 


Previous Theoretical Results and Experimental Goals 
This section summarizes the previous theoretical results in 
relation to the goals of current experimental work. 

When both the number of structural modes and the number 
of acoustic modes are large then, from Peretti and Dowell 
(1992a) the nondimensionalized mean square pressure locally 
(i.c.* as a function of position in the cavity (can be expressed 


as: 
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From, (I), the quantities which are needed to predict the local 
piean-square pressure are: the flexible wall area (Aj), the vol- 
ume of the cavity ( V ), the acoustic modal function, (F r (.v,j\z), 
cosines for the rectangular cavity), as well as quantities which 
are evaluated at the center frequency of the band. The acoustic 
generalized mass (M C A ), the center frequency (« c ), the acoustic 
modal damping (f/), and the power spectra of the wall ac- 
celeration ($ w ) are all evaluated at the center frequency only 
and are not determined for each individual mode. F c stands 
for the modal average of the acoustic modal function. (Z f >is 
the ratio of the volume average of F c 2 to the average of F c 
over the flexible wall area. The following assumptions are 
implicit in this equation: the power spectra of the wall response 
is uncorrelated in space, which is a consequence of having a 
large number of structural modes; and the power spectrum of 
the wall response is slowly varying with respect to frequency 
relative to the rapidly varying acoustic cavity transfer function. 
Also, it is assumed that the generalized mass, the frequency 
in the band of interest, and the damping ratio vary slowly, 
and can be approximately by their values at the at the center 
frequency of the band. A complete derivation of Eq. (1) is 
given in the Appendix. 

Increasing the bandwidth, and therefore increasing the num- 
ber of modes at a given center frequency, should improve the 
agreement between t lie experimental results and the asymptotic 
modal analysis results. However, beyond a certain bandwidth 
size, the assumptions that frequency, damping, etc. are con- 
stant are violated, and this tends to negate the positive effect 
on analytical/experimental correlation which is due to the in- 
creased number of modes. One goal of the experiment was to 
investigate the bandwidth effect. 

When there is a large number of acoustic modes in the cavity, 
i.e., when the index r in Eq. (1) goes to infinity, the mean 
square pressure levels in the interior approach a uniform value. 
The only term in the AMA expression for mean square pressure 
which contains spatial variation is inside the summation over 
acoustic modes. However, for very large r, the summation of 
F f l equals a constant, which is the same for all x, y % and z, 
unless .v, y , or z are close (within an acoustic wavelength) to 
the cavity boundaries. The condition, x, and/or y , and/or z 
equal to 0, L T , L yy or L zt corresponds to a wall, an edge or a 
corner location. , 

For a rectangular enclosure, substituting Fj = cosine _ (r t 7r.v/ 
L x ) * cosine 2 (r,iry/Z^) * cosinc\r z vz/LJ into Eq. (1), and setting 
either x, y t or z equal to 0 or L x , L yy or L Zi for a wall, the 
term LF^/AN* for large r becomes 1/4. In the interior region, 
where neither x, y , nor z are equal to 0 or the cavity dimensions 
in their respective directions, this ratio is 1/8. Therefore, the 
ratio of wall mean square pressure to the uniform interior 
pressure is 2 to 1. Similarly, the ratio of mean square pressure 
in a corner to the uniform level is 8 to 1, and for an edge it 
is 4 to 1. These intensification zones at the boundaries are due 
to the spatial correlation of the modes there (enforced by the 
rigid wall boundary conditions) as opposed to the random 
distribution (spatially and temporally) in the interior. Vaii- 
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dation of the uniformity of the interior sound level as well as 
the intensification levels at the walls, edges and corners was 
another goal of the experiment. 

The intensification zones exhibit boundary layer behavior, 
and are relatively flat in the interior region. Peretti and Dowell 
(1992b) studied the transition from the peak boundary levels 
to the uniform interior level for the rectangular acoustic cavity 
previously considered. They define the nondimensiona! pa- 
rameter, k<d 9 as the center frequency wavenumber multiplied 
by the distance into the cavity. The mean square pressure as 
a function of k<d is independent of cavity dimensions and is 
only dependent on the bandwidth to center frequency ratio 
and direction from the cavity boundary. The thickness of the 
intensification zone scales with the center frequency wave- 
number and is such that for different bandwidth to center 
frequency ratios the intensification thickness is relatively con- 
stant. A further goal of the experiment was to measure the 
sound pressure levels as they transition from the peak levels 
at the boundaries to the nearly uniform levels inside the cavity. 

Experimental Design and Description of Apparatus 

The experiment was designed to match as closely as possible 
the theoretical assumptions. These assumptions include: a rigid 
rectangular cavity with only a flexible portion of of one w'all 
vibrating, a large (infinite) number of structural modes re- 
sponding on the vibrating wall portion, and a large number 
of acoustic modes containing all frequencies in the band (“white 
noise”)- In the theoretical work, the width, height, and length 
dimensions were in the proportions 2 to 3 to 7. These dimen- 
sions were chosen so that they would not be integer multiples 
of each other. 

A rectangular box was constructed almost entirely of 1 in. 
thick Plexiglass. The interior dimensions of the box were 408 
mm X 610 mm X 1520 mm. One or the 408 mm by 610 mm 
walls was constructed of Bakelite and was removable. This 
feature allows the flexible plate location and dimensions to be 
varied in subsequent experiments. The portion of the Bakelite 
wall that was rigid was 25,4 mm (1 in.) thick, and the thickness 
of the portion which was flexible varied according to its other 
dimensions. The requirement of having a large number of 
structural modes responding on the flexible Bakelite plate dic- 
tated the thickness to area ratio of the plate; the smaller the 
plate area, the thinner the plate needed to be. For the exper- 
imental data reported here, the flexible area was positioned in 
the lower left corner of the wall, its area was 622.2 square cm, 
which corresponds to 1/4 of the entire wall area for that face, 
and its thickness was i.5875 mm (1/16 in.). 

The flexible portion was driven by white noise which was 


supplied by a B&K Noise Generator connected to an amplifier 
and then to a 2 lb. shaker. The sting of the shaker was attached 
to the flexible portion of the wall by wax at a position which 
was off-center, so as to excite both symmetric and antisym- 
metric modes. Previous work by Kubota and Dowell (1986) 
has shown that when a plate is driven in this manner, a large 
number of (structural) modes respond, and the response char- 
acteristics are nearly spatially uniform; exceptions occur at the 
structural boundaries and at the point of excitation. For this 
type of behavior, at least fifteen structural modes must be 
present in the bandwidth. Figure I shows a plot of the number 
of structural modes for the flexible plate used in this study as 
well as the number of acoustic modes in the cavity versus 
frequency for various frequency bandwidths divided by center 
frequency, i.e, .005, .2329, .500. 

An accelerometer was used to measure the acceleration of 
the vibrating plate. The accelerometer was carefully positioned 
so as not to be near a boundary or the force application point. 

Microphone probe tubes were used to take sound pressure 
level measurements inside the cavity and they were approxi- 
mately 20 cm (8 in.) long. Holes of approximately 4 mm in 
diameter were drilled at various locations on the 5 rigid walls 
to accommodate the microphone probe tubes. In order to keep 
the cavity “rigid,” the holes where plugged vyhen not in use. 
The end of each plug was aligned with the interior of the cavity. 
A schematic drawing of the cavity is shown in Fig. 2. 

The microphone probe tubespossessed standing waves which 
were damped out using open cell foam strips, which were 
inserted into each tube. Figure 3 shows a diagram of a micro- 
phone probe tube. The effect of the damping was accounted 
for in the final pressure measurements by calibrating accord- 
ingly. The calibration was done by comparing the outputs of 
two microphones: a half-inch microphone without a probe 
tube, and a half-inch microphone with a probe lube. The two 
microphones were placed in front of a loudspeaker which was 
driven with white noise, such at that the mouth of the probe 
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Fig. 3 Schematic of microphone probe tube 
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Fig. 4 Transfer function of the response of microphone probe tube 
without damping to a bare microphone, and transfer function of the 
response of microphone probe tube with damping to a bare microphone 


tube was flush with the face of the bare microphone. The 
transfer function between the two microphone outputs was 
taken on a spectrum analyzer and stored on a computer for 
future use in calibrating the data. Figure 4 shows transfer 
functions* with and without damping. The added damping has 
the effect of removing the standing wave peaks without a large 
reduction in signal to noise ratio. 

The ratio of the output (in volts) from the two microphones 
(without probe tubes) subject to a fixed pressure input was 
also needed for the calibration. This ratio was assumed to be 
independent of frequency in the frequency range of interest, 
and therefore, could be obtained from two pistonphone meas- 
urements. Since the system is linear, the raw data could be 
corrected for the damping effect in the probe tube by using 
the ratio of the two microphones without probe tubes, and the 
transfer function (which varied with frequency) for a micro- 
phone with a probe tube to one without. 

During the experiment, two microphone/probe setups (both 
calibrated as described earlier) were used to measure sound 
pressure levels at two interior points simultaneously while the 
accelerometer measured the motion of the wall. These three 
data samples were input to a spectrum analyzer (4-channel 
Scientific Atlanta SD-380), from which power spectra could 
be plotted on a compatible Hewlett-Packard plotter directly, 
or the raw data could be stored on the VAX computer via an 
IEEE connection, and later post-processed. The plotter option 
was convenient for checking measurements in the early stages 
of the experiment. However, the bulk of the experimental data 
was stored on a computer as raw data (i.e., volts) and later 




Fig. 5 Experimental set up schematic 


post-processed, since the goal was to compare w’ith the earlier 
analytical study in which nondimensional groupings were used. 

Results and Discussion 

Experimental Procedure and Data Process- 
ing. Microphones with probe tubes were used to measure 
sound pressure levels at different locations inside the rectan- 
gular acoustic cavity while an accelerometer measured the 
shaker input to the wall. For each fixed microphone probe 
tube location, three independent readings were taken for both 
the microphones and the accelerometer. The microphone and 
accelerometer data were read into a spectrum analyzer and 
converted to power spectra densities. The power spectra ranged 
from 0 to 5000 Hz and consisted of 400 data points. Therefore, 
the narrow bandwidths for the raw data were 12.5 Hz. The 
power spectra data in volts were stored on a computer for later 
post-processing. 

The first step in the post-processing was to read in the three 
independent measurements for a fixed spatial coordinate inside 
the cavity and average them. The measurements were checked 
by plotting the three samples and their average. These data 
were then converted from volts into mean square pressures in 
Newtons per meter squared by using the earlier microphone 
calibration data and the conversion factors from the piston- 
phone measurements . The second post-processing step was to 
take the average of the three narrow band power spectral 
densities (12.5 Hz spacing), and create new spectra with dif- 
ferent bandwidth to center frequency ratios. For this, a com- 
puter program was written which stepped through center 
frequency, computed the bandwidth for a desired bandwidth 
to center frequency ratio, and summed the squares of the data 
in that desired band. The outputs of this computer program 
were power spectral densities (in dB) in continous bands, typ- 
ically from 200 Hz to 4000 Hz w'ith spacing every 20 Hz, for 
a particular bandwidth to center frequency ratio (/«//<•). Typical 
plots for a microphone and an accelerometer are shown in Fig. 
6. 

The third and final step in post-processing introduced the 
probe location into the data. Here, a series of power spectral 
densities for a given f b /f c ratio were read into a computer 
program, along with their associated position (x, v, and z 
coordinates) in the cavity. The output of this program was 
mean square pressure level versus k t d, where k^d is the wave- 
number (fc c ) for a particular center frequency times the distance 
(d) into the cavity with respect to a reference point. Upon 
completion of the final step, the data have been transformed 
from power spectral densities in volts to mean-square pressures 
as a function of a nondimensional parameter k,d for a given 
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Fig. 6 Typical accelerometer and microphone data for psd's with 1/3 
octave band 
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Fig. 7 Scatter of the psd data at 6 different locations away from the 
boundaries of the cavity. Numbers in parenthesis indicate the number 
of modes in a 1/3 octave band at that center frequency 
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center frequency and given bandwidth to center frequency 
ratio. 

Uniform Interior. When there is a large number of spatially 
uncorrelated acoustic modes in the cavity, the theory predicts 
that the interior sound pressure levels in a band will be uni- 
formly spatially distributed. To test this hypothesis, micro- 
phone probe tubes were inserted into the cavity to their 
maximum extension (approx. 20 cm) at various locations on 
the five rigid walls. The power spectral densities (in dB re. 2 
X 10" 5 N/M 2 ) were plotted versus frequency (Hz) for a 1/3 
octave band. The scatter of this data is shown in Fig. 7, where 
the deviation from the mean is plotted versus frequency. Note 
that above 2000 Hz, the data contain little scatter, while below 
2000 Hz the scatter in the data is as much as ±10 dB. This 
suggests that the region above 2000 Hz will be in the AMA 
limit, that is, the theoretical conditions are met regarding the 
large number of acoustic modes. For this geometrical config- 
uration, a frequency of 2000 Hz corresponds to 26 plate (struc- 
tural) modes and 266 cavity (acoustic modes in a 1/3 octave 
band. This compares to 16 plate modes and 64 cavity modes 
at 1000 Hz. (Figure 1 shows the growth in the number of modes 
versus frequency for various bandwidth to center Irequency 
ratios.) 


Intensification Zones Near The Cavity Bounda- 
ries. Theoretically, response peaks occur at the boundaries 
of the cavity, due to the imposition of the boundary conditions. 

These peaks are 8, 4, and 2 times higher than the spatially 
uniform values that exist further into the cavity; the level of 
intensification depends upon whether the boundary point is 
on a corner, edge or a wall, respectively. Analytically, these 
intensification zones have been studied by Peretti and Dowell 
(1992b) by plotting the mean-square pressure as a function of 
nondimensional distance into the cavity (kcd). The two main 
points for comparison between theory and experiment for the 
intensification zone study were the ratio of the levels at the 
boundaries to those in the uniform region and the thickness 
of the intensification zone (as determined by the distance to 
the first minimum). 

Figures 8, 9, and 10 show intensification curves coming from 
a corner, an edge, and a wall into the interior. In the cases of 
the corner and the edge, the trajectory was at an oblique angle, 
whereas for the wall, the trajectory was perpendicular to the 
wall. The three cases shown here are for a 1/3 octave band- 
width. Several center frequencies are plotted for the experi- 
mental data; the theoretical data is only dependent upon the 
ratio of bandwidth to center frequency which is .2^29 tor a - - 

1/3 octave band. The ordinate is the ratio of mean-square 
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Fig. 10 Intensification curves for the wall to interior case 


pressure at the x , y , z point to mean-square pressure of the 
interior region. The abscissa is the nondimensiona! distance 
into the cavity. 

The theoretical curve in Fig. 8, for the case where k,d starts 
in the corner (0, 0, 0) and proceeds into the cavity at an angle 
which is approximately 33 deg in one direction and 56° in the 
other, starts with a ratio of 8 in the corner, and then becomes 
flat with a level of 1 for large values of k<d. Also, the first 
minimum occurs around k<d = t, which corresponds to d * 
X<V2, where \ c is the center frequency wavelength. The exper- 
imental curves which best fit the theoretical curve are for fre- 
quencies above 2000 Hz. This was expected since 20000 Hz 
seemed to be the lower bound on the AMA regime from the 
earlier results evaluating the uniformity of the interior sound 
levels. Therefore, at 2000 Hz, there are not enough modes for 
the theory to hold, although the agreement is still fairly good. 
The experimental data show better agreement with theory as 
the center frequency and, therefore, the number of inodes 
increases, up to a point, and then it is slightly worse. At 4000 


Hz, the ratio of mean-square pressure in a corner to that in 
the interior is almost JO to I. The deterioration of the theory 
at higher center frequency is due to the fact that the bandwidth 
is proportional to the center frequency. Therefore the band- 
width gets wider as the center frequency. Therefore the band- 
width gets wider as the center frequency increases. In large 
bandwidths, the assumptions that the frequency, power spectra 
of the wall acceleration, damping, etc. are essentially constant 
may break down. 

Similar results are shown in Figs. 9 and 10. In Fig. 9, the 
theoretical curve has a value of 4 at k <d equal 0(at the edge), 
reaches its first minimum around 3 oscillates slightly and 
asymptotes to 1 for large values of k t d. Here, k<d is the non- 
dimensional distance from the edge along a line which makes 
an angle of 33 deg. with the wall. Experimental results agree 
well for center frequencies above 3000 Hz. At 3000 Hz, the 
value of mean-square pressure at the edge is more than 5 times 
higher than the interior. 

The mean-square pressure ratio in the intensification zone 
near a wall (away from its corners or edges) is shown in Fig. 
10. Theoretically, the ratio is 2 to 1 at the wall. Again, agree- 
ment is good, although the experimental levels are higher than 
the predicted levels at the wall, for the highest center frequency 
values. 

In all three cases, the experimental data show asymptotic 
behavior, that is, the mean-square pressure levels approach a 
constant for large k<d. In almost all instances, the theory pre- 
dicts the thickness of the intensification zone well. 

Effect Of Bandwidth. The previous data suggest that the 
theory holds when the center frequency is large such that there 
is a large number of acoustic modes in the cavity. However, 
it also suggests that when the bandwidth is too large the theory 
breaks down. In Fig. 1 1 , the center frequency is held constant 
at 3000 Hz, and the bandwidth to center frequency ratio is 
varied from 5 percent to 50 percent. The experimental data 
which are used here are for the wall intensification zone. Six 
experimental curves are shown. A corresponding theoretical 
curve for each experimetnal curve is not shown in order not 
to crowd the plot. Although there arc six different theoretical 
curves to compare to, all the theoretical curves have a value 
of two at the wall, and asymptote to one away from the wall. 
Therefore, qualitatively this data can be judged on how well 
it agrees at the wall and asymptotically. The worst cases are 
the data corresponding to the 5 percent and the 50 percent J h 
/•ratios, but for different reasons. For small bandwidths, there 
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Fig. 12 Damping coefficient, f/, predicted from reverberation time 
measurements 
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Fig. 13 Comparison of experimental and theoretically predicted values 
for spatially-averaged mean-square pressure 


may not be enough modes to satisfy the condition that the 
number of acoustic modes be very large. However, for large 
bandwidihs, the condition that certain variables be considered 
constant in the band is violated. 


Prediction of Interior Levels using AMA vs. Experimental 
Results. The spatially-averaged mean-square pressure of the 
interior is described by AMA theory (see the Appendix) as: 


<p 2 > x A N** w 
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The quantities A s (the flexible wall area) and V (the volume 
of the interior) are known from the geometry of the rectangular 
acoustic cavity. The number of modes AN* in the desired 
bandwidth can be estimated using the standard formula from 
architectural acoustics (in any of the textbooks referenced). 
The term <Z C 2 > is a ratio of the acoustic modal functions 
evaluated at the center frequency spatially averaged over the 
entire volume to that spatially averaged over the flexible wall 
area. For this example, the acoustic modal functions are prod- 
ucts of cosines. The damping ratio evaluated at the center 
frequency, f/ can be estimated from independent measure- 
ments of the reverberation time in the acoustic cavity, which 
were performed as part of the overall experimental effort. 

The reverberation time was measured using a B&K sound 
level meter with a reverbation module attached, (B&K Sound 
Level Meter Type 2231 Plus Reverberation Processor Module 
BZ 7104) which sends 1/3 octave band bursts and then meas- 
ures decay as a function of time. It calculates and displays the 
reverberation time in each 1/3 octave band. The damping ratio 
was predicted from reverberation time measurements through 
the relationship described by Dowell (1978) and are plotted in 
Fig. 12. The remaining unknown parameter in Eq. (2) is the 
power spectra of the wall acceleration, <f>*. This information 
can be obtained from the accelerometer measurements shown, 
for example, in Fig. 6, which were taken concurrently with the 
microphone readings for sound pressure levels. The acceler- 
ometer was moved around on the flexible portion of the wall, 
and the response was found to be nearly spatially uniform 
above 2000 Hz, i.e., the flexible plate, itself, was in the AMA 
limit for structural vibration above 2000 Hz (approximately 
26 modes in "a 1/3 octave band). 

A comparison between experimentally obtained results and 
theoretically predicted values for spatially averaged mean 


square pressures are shown in Fig. 13, The solid line corre- 
sponds to data which were measured experimentally with mi- 
crophones and then averaged. The spatially-averaged sound 
pressurejevels which are predicted using AMA, with the input 
to the AMA formula being the accelerometer measurement 
from the flexible wall and the predicted damping ratio from 
reverberation time measurements, arc represented by discrete 
points. An additional curve is also shown in the figure, it 
corresponds to the AMA predicted values for the mean-square 
pressure assuming a constant damping ratio of 3 percent. The 
agreement shown is quite good, particularly above 2000 Hz. 
where the flexible plate is in the AMA limit. Below this value, 
one accelerometer reading is not sufficient to represent the 
power spectra of the wall acceleration. It is believed that better 
agreement could be obtained by using a spatial average of the 
output from several accelerometers on the flexible wall for 
<f>*. However, below a certain value of center frequency, the 
assumptions for AMA, both in the acoustics as well as the 
structural dynamics, break down and good agreement can never 
be reached. 


Conclusions 

When there is a large number of structure modes responding 
on the vibrating portion of the wall and when there is a large 
number of acoustic modes in the cavity, the AMA theory 
predicts that the sound pressure levels in the interior are uni- 
form. This was verified by the experiment performed here, 
which showed that for this cavity configuration, an AMA 
regime exists above 2000 Hz for power spectral densities taken 
in 1/3 octave bands. The AMA theory also predicts that in- 
tensification occurs at the cavity boundaries, and that the re- 
spective ratios relative to the interior level are 8, 4, and 3 for 
corner, edge, and wall intensification. The theory and exper- 
iment showed reasonable agreement with respect to the levels 
of intensification at the boundaries as well as the thickness of 
the intensification zone. The combination of bandwidth and 
center frequency must be carefully chosen to give optimal 
results from the AMA theory. If the bandwidth is too narrow, 
or the center frequency too low there may not be enough modes 
for the assumption of infinite number of acoustic modes to 
be valid. Whereas, if the bandwidth is too wide, other AMA 
assumptions break down, such as that the modal properties 
are slowly varying over the bandwidth. 
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APPENDIX 

The following is a derivation of Eq. (1) in the text. Equation 
(2) is a spatial average of Eq. (1). A Similar derivation can be 
found in Kubota, Dionne and Dowell (1988). 

The structural wall modal equation of motion is 

MmlQm + 2tmumq m + u m 2 q„,] = Qj ( 1 ) 

where the modal expansion for the wall deflection is 

w=^q „,(!)■* „(x.y) ( 2 ) 

m 

and the structural generalized mass is 

| nipt/ m \x,y)dxdy ( 3 ) 

and the generalized force due to a given external pressure is 
Qj* j j p E V m (x,y)dxdy (4) 

The acoustic cavity modal equation is 

Pr + Utr Pr + r~Pr = Qr (5) . 

where the modal expansion for the acoustic cavity pressure is 


2 Pr{t)F r(-YJ',Z) 

p=p *° 


( 6 ) 


(7) 


and the acoustic generalized mass is 

*-k!U F,\xy,z)dxdydz 
and the generalized acceleration due to the structural wall is 
Q ? ’■•p j | j wF r (x,y,z)dxdy. (8) 

Define /, a nondimensional cavity pressure, 

P 


At,xj,z)* 


Po^o 


From (6) and the defintion of auto-power spectrum, the auto- 
power spectrum of / is: 


^ / ^ FVCv,>',z) F 5 (xJ,Z) A 

*/W,z)=2j Zj ~w' mT 


where the cross-spectra are defined as: 


<f> 






(ry-'aV 


( 9 ) 


(10) 


and the cross-correlations of the modal generalized pressure 
coordinates are: 

i r r 

— P,(t)PAt + T)dT. (II) 

7— OB Z 1 J __ f 

Similarly from (8) the cross-power spectra Q / and Q s " are: 
<t V , c J H (“> = p jjjj F r (x,y,z)F s (x‘ ,y’ ,z‘)‘t>,(u;x,yy ,y’ ) 

dxdydx* dy' . (12) 

The transfer function for the acoustic modal Eq. (5) is: 

1 


^rM- A j. 

— 0) + 2 lifter W + (jJ r 


( 13 ) 


From standard random response theory and Eq. (5), the re- 
lationship between $> P/fis and $ r *Q s w is: 

*^,(o>) = - u)* Q r Q »iw). (14) 

From (9). (12). and (14): 

.1 x-i vi FXx,y,z) F s (x,y,z) r , A , ,, tA , 

4>/ u pr j-.z) = p 2 j Zj —jT* Hr )H « ( ~ w) 

r s r * 

( 15 ) 


'HI | F r (x,y,z)F s (x‘ ,y' ,z * )4> h .(cj \x,y,x ' ,y * )dxdydx ' dy * 

This is the basic expression for the power spectra of the 
cavity pressure in terms of the power spectra of the wall ac- 
celeration. The double summation over r and s can be reduced 
to a single summation by ignoring the off-diagona! coupling. 
This is called ‘‘joint acceptance” theory, and it entails omitting 
all terms in the double sum except those for which r = 5, This 
theory is applicable because - w) « ! ! 2 and H s 

{w)-H r (-o>)«\H s (o>)\ 2 for typical small damping values. 
Equation (15) is then simplified to: 




l y 1 

T 3 ^ (A !?? 


l// f 4 (aj)l 2 


( 15 *) 


•JJJ | F r {x,y\z)F f (x' y ,c*)^ (u.';.v , v 

Since, 


)dxdydx* dy' 


2 V m {xj)ty n {x\y')w~H n XufH m (uWH n ( - < 


2 *Jx l <J , ,<*n(x/j/)*r 0 (u) 

U 

then, for a large number of excited structural modes ( AM — co), 

, , 0 ifx*x* 

Q^xjpc j )& 

constant if(jc = jc ,y-y 

There f ore , $ w (u> x* y) = Aft k .(u>)5(x - x* )5(y - > * ). 

Assuming the power spectra is smoothly varying with respect 
to the transfer function (“white noise assumption’*), then 
$ w (u)s<l >n(w f ). Equation (15*) then becomes: 


The integral over the flexible area can be expressed as a 
spatial average over the flexible area limes the area, Fj. The 
acoustic modal function, F n can be written in terms of spatially 
independent functions, X r (x)*Y r (y)*Z r (z). Assuming that the 
Flexible wall area is in the z = z 0 plane, 
F r {x % y>Zo) = XM Y r (y)Z r (z 0 ). Mathematically, 


{{/*«. 


z a )dxdy 


<^xy;xz;) _ <^> 

Volume 

<Zf) <Zl) ' 


■I L 


)dxdy ( 1 6 ) 


Integrating over frequency to obtain mean-square pressure 
gives: 

72 P 2 _ y Fl(x,y } z) 

*~W> 2 " 4 K 2 


* F;(x,y f Zo)dxdy_ (17) 


This is the equation for mean-square pressure in the cavity 
when there is a large number of structural modes and a mod- 
erate number of acoustic modes. Assuming A/?, w?, and 
do not vary rapidly in the bandwidth, and can be replaced by 
their values at the center frequency of the band (subscripted 
by “c”), Eq. (17) becomes: 

V F:(r v r) 


(^7? 4 P ( A/^V) y (F) ^ 


When there is a large number of acoustic modes, (r— co) 
and when the acoustic wavelengths are small compared to the 
size of the flexible area, then </^> and <Z 2 > are independent 
of modal number, r, and can be taken outside the summation. 
If these conditions do not apply, then </*?) and <Z;> should 
be modally averaged over the frequency band. Subscripting F 
and Z by “c” to denote that they are independent of r, and 
taking them outside the summation, Eq. (19) becomes: 

0yV) 2 ~4 V 1 <Z 2 > 2 ' ( *** z) ( > 

Equation (20) is identical to Eq. (I) in the text. Equation (2) 
in the main text is simply a spatial average of this equation. 
Using the following relations, 

= and J] </>(.vj',c)> S AA^.</^> 

r 

in the above equation leads to the spatially averaged result: 
<P 2 > . T (aA 1 AN^ 


j F i r {x,y,z 0 )dxdy (18) 


ifiaCoT 4 \ V (o-v) J ^<Zr) 


which is the same as Eq. (2) in the text. 
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CHAPTER 4 


Asymptotic Modal Analysis for Structural-Acoustic Systems with Finite 
Structural and Infinite Acoustic Modes 


Introduction 

When part of an enclosure vibrates, it often excites a large number of acoustic 
modes in its interior. Such is the case when part of an automobile vibrates, creating 
unpleasant acoustic responses in the interior, or when the sidewall of an airplane fuselage 
vibrates, transmitting annoying sound pressure levels to the interior. The results presented 
here are directly applicable to these vehicular noise and other interior noise problems. 

Previously, Asymptotic Modal Analysis has been used to study a similar class of 
problems in structural-acoustics [Ref. 4, 11, 22, 23]. In the previous works, it was 
always assumed that the number of structural modes was large, and the number of acoustic 
modes was either large or small. Here, it is assumed that the number of acoustic modes is 
large (approaches infinity), and the result will be derived for the case where the number of 
structural modes is finite. 

This result, in the limit where the number of structural modes approaches infinity, 
will be compared to the result derived previously for an infinite number of both structural 
and acoustic modes. In the previous result, it was initially assumed that the number of 
structural modes was infinite, and then the number of acoustic modes was allowed to go to 
infinity. 

Asymptotic Modal Analysis (AMA) is a recent method, which is systematically 

derived from classical modal analysis. AMA was originally developed by Dowell [1] as a 

way to demonstrate the relationship between statistical energy analysis and classical modal 

analysis. The book by Richard Lyon [5] is the classic reference on statistical energy 
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analysis. The AMA method has many advantages. Since it is systematically derived from 
CMA, the generality of the final result can be modified by changing assumptions which are 
made in the derivation. For instance, AMA can predict spatially averaged sound pressure 
levels in the interior of an acoustic cavity, and it also has the capability to predict local 
sound pressure levels. The AMA method requires less information about the system than a 
traditional CMA approach. Parameters, such as modal damping ratio, frequency, 
generalized mass, etc. for each individual mode are needed for CMA, whereas in AMA, 
only the values at the center frequency are required. 

The AMA method is useful for solving linear dynamical systems where there are a 
large number of modes and has been proven accurate for both structural dynamics and 
acoustics problems. The method has been verified experimentally by Kubota [3] for a 
rectangular plate excited by point forces as well as a vibrating plate with a concentrated 
mass or masses [24, 25]. Experimental verification of the AMA method applied to a 
rectangular acoustic cavity in which part of one wall vibrated was reported by Peretti & 
Dowell [22], In all of the previous AMA analyses it was assumed that the number of 
structural modes was large. Conversely, for the structural-acoustic system studied here, a 
small number of structural modes will be assumed initially. The underlying assumption in 
this analysis is that the number of acoustic modes approaches infinity. 
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Theoretical Development 


Derivation of the Result for a Finite Number of Structural Mod es, and an Infinite Number 
of Acoustic Modes 

The expression for mean-square pressure as a function of wall acceleration for a 
rigid wall acoustic cavity with a flexible vibrating portion is derived below. This 
derivation is similar to the derivation found in Ref. 4. However, in Ref. 4 the initial 
assumption is that the number of structural modes is infinite . Here, the underlying 
assumption is that the number of acoustic modes is infinite . The derivations are identical 
up to Equation (4.14), but all of the steps are included here for completeness. 

The structural wall modal equation of motion is 

M Jq m +2CmCOmq m +CO m 2 q m ] = Qm E (4 ’ 1 ) 

where the modal expansion for the wall deflection is 


w = £q m (t) YJx.y) 

m 

and the structural generalized mass is 

-II 


M n 


m p' F m fr.y) 




and the generalized force due to a given external pressure is 


Qm 1 


■ 11 ' 


p Ev F m (x,y) dxdy 


(4.2) 


(4.3) 


(4.4) 


The acoustic cavity modal equation is 

A A. A vu (4 5) 

p r +2£ r co r p r +co r p r = Q r 

where the modal expansion for the acoustic cavity pressure is 

lx' PrtoFiOtiy.z) 

P = Po c o2, -T (4.6) 

r M r 
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and the acoustic generalized mass is 

Mf = y//J F r 2 (x,y,z) dxdydz 

and the generalized acceleration due to the structural wall is 

wFj(x,y,z) dxdy . 

V A f 

Define f, a non-dimensional cavity pressure. 


f(t,x,y,z) s — 


(4.7) 


(4.8) 


Po C o 

From (4.6) and the definition of auto-power spectrum, the auto-power spectrum of f is: 

= (4.9) 

where the cross-spectra are defined as: 

Rrps ( t ) etot dx 


^PtPs - 


■if. 


(4.10) 


and the cross-correlations of the modal generalized pressure coordinates are: 

r ivps = lim Pr(t)P s (t+t)dT . 


£« 


(4.11) 


T-*~ 2Tj 

Similarly from (4.8) the cross-power spectra of Qr w and Qs w are: 

Oq^co) = ^JJJjF r (x,y,z)F s (x*,y*,z*)Ow(co;x,y,x*,y*) dxdydx*dy*. (4. 12) 

The transfer function for the acoustic modal equation (4.5) is: 


H r A (co)=- 


1 


2 A A 

-co +2i£,co r co+co r 


(4.13) 


From standard random response theory and equation (4.5), the relationship between Op r p s 


O ftR (co) = H A (co)H A (-co) Oqw^(co) . 


(4.14) 


From (4.9), (4.12), and (4.14): 


<D f (co;x,y,z) 


_l_yy Fr(x ’ y ’ z) M H A (co)rf(-co) 

W 2 T 7 H A M A 


(4.15) 


F r (x,y,z)F s (x*,y*,z*)O w ((o;x,y,x*,y*) dxdydx*dy !| 


This is the basic expression for the power spectra of the cavity pressure in terms of 
the power spectra of the wall acceleration. The double summation over r and s can be 
reduced to a single summation by ignoring the off-diagonal coupling. This is called joint 
acceptance" theory, and it entails omitting all terms in the double sum except those for 
which r=s. This theory is applicable because H r (co)H s (-co)« IH r (co)l 2 and H s (co)H r (- 
co)« IH s (co)l 2 for typical damping values. Equation (4.15) is then simplified to: 


® f (o*x.y.z) = M : 

V 2 r (M A r 


•Ik y,z)F I (x !k ,y*,z*)Ow(0);x,y,x*,y*) dxdydx*dy< 


(4.15*) 


Further, w = Xq m (t)Y m (x,y) , from (2) 

m 

Rw= lim w,(x,y,t)w s (x*,y*,t+T) dt 

= lim -l-[ Iq m (t) v F rn (x,y)Idn(t + x)T n (x*,y*)dt 
T-»co 2TJ t m n 


or, R^(x,x,y,x*,y*) = £ £ ^ • 

m n 

Taking a Fourier transform (t — > co , R — > O) 

®rfco,x,y,x*,y*) = SSW'P.fxWii.W (4.16) 

m n 
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Again, ignoring the off-diagonal coupling by "joint acceptance" arguments, and noting that 
the point pairs x, y and x*, y* correspond to the same points and are therefore, 
superfluous: 


|H^m)| 2 -[JJ F 1 (x,y,zJ'P4x,y)dxdy] 2 -4> w ico) (4.17) 
V 2 r (mA) 2 L Af J 


Equation (4.17) is for one structural mode (m). The total mean- square pressure can be 
computed by summing the individual contribution from each mode. The mean square 
pressure response is found by: 


(poC^f 


c=0;x,y,z) = j 


n 2 

f 2 = — 2 — = Rj(t=0;x,y,z) = I Oj(co;x,y,z) dco 


Integrating equation (4.17) with respect to co, and recalling that the white noise assumption 


implies that d>- • is relatively constant in the bandwidth, yields: 

r qrrQt 


p 2 - 1 it y F, 2 (x.y,z) i 
(p 4 = v i4 T (ma ) 2 UK 


■ [JjAf F i( x >y> z o)^m(x,y) dxdyj- O^co) (4. 1 8) 


If the flexible area is small compared to an acoustic wavelength, F r (x,y,Zo) can be 
replaced by its value at the centroid, F r (x 0 ,yo> z o)- The implication and justification for this 
assumption will be discussed later. The integral over '?m(x,y) is equivalent to the spatial 
average of x F m (x,y) times the area of the flexible wall (Af). Equation (4.18) can then be 

rewritten: 


P 2 -x(Af) 2 V F ^ (x,y ' 2) •{'F„(x,y)K f * (4.19) 

(poC 2 ) 2 ” 4 V ^ (MA ) 2 (coAfcA 

Since, there are an infinite number of acoustic modes (r-*») the summation over r 
can be approximated by: 

y = LMMm I . o w (co c ) (4 . 20 ) 
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This assumes that the generalized mass (M r A ), damping (Cr A )> frequency (03 r A ), and power 
spectra (O(C0r A )) are slowly varying with respect to the rapidly varying transfer function 
(H r A ) and can therefore be replaced by their values at the center frequency. The summation 
over the product of the acoustic modal functions is then replaced by the product of their 
spatial averages times the number of acoustic modes in the band. 

Using these approximations in equation (4.19), the result for the mean-square 
pressure when there are a large number of acoustic modes, but a finite number of structural 
modes, is represented by equation (4.21). 


P 2 


(poCo 2 ) 2 


-S-fM 2 faf (x y) )af (Fc)°AN . Qq; o’ (co c ) 

■ 4 lv) 1 m( ,y)/Xf (f^Cc 


(4.21) 


This is the result for the mean square pressure when one structural mode excites a large 
number of acoustic modes in the bandwidth. The total mean-square pressure in the band is 
the summation of the contributions from each structural mode in the band. 

This expression can also be written as a function of the power spectra of the 
individual point forces rather than the power spectra of the wall acceleration in generalized 
coordinates. That result and its derivation are presented in the appendix. 


Result for an Infinite Number of Structural Modes and an In finite Number of Acoustic. 
Modes 

Equation (4.21) is the result for the mean-square pressure when there are an infinite 
number of acoustic modes excited by one structural mode. It is instructive to take the limit 
of this expression summed over all structural modes, m, as m goes to infinity. This will 
produce a result which is valid for a case where there are both a large number of acoustic 
and a large number of structural modes in the band and can be compared to the previously 
derived result [Ref. 4] where the limits were taken in the opposite order. 
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In the AMA structural mode limit, ('I'm) 2 — > (p) 2 , and ('I'm) -> ('F 2 ) an( ^ 

_ (» 2 ) 

— / o \ , which follows from (4.16). Using the above relations and 

yFc / AM Aco 

summing (4.21) over m, yields: 

p 2 _ tt / AfY 2 an a 1 ('Pcl't.y )? 

(p 0 c If 4 V AO) (tOcpCo ('P 2 (x,y)> 

where <Zc 2 > 0 is the ratio of the spatial average of F c 2 over the volume to F c 2 at the centroid 
of Af . This is the result for the mean-square pressure when it is first assumed that an 
infinite number of acoustic modes responds, and then that an infinite number of structural 
modes also responds. 

Previously, Kubota, Dionne and Dowell [Ref.4] derived a similar result for an 
infinite number of both structural and acoustic modes. Their derivation first assumed that 
an infinite number of structural modes was responding and then took the limit as the 
number of acoustic modes approached infinity. From Ref.4, that result is: 

(p 2 ) B s. /Afj 2 AN (w^acq (4-23) 

(PoCo 2 ) 2 _ 4 V A “ A (^FCc A (4 

The two results are seen to be very similar. From dimensional considerations this 
is hardly surprising. What is significant are the differences. A discussion of the 
differences between these two asymptotic results follows. 




(4.22) 
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Discussion of Results 


Equations (4.22) and (4.23) show that the result for an infinite number of acoustic 
modes and an infinite number of structural modes is not independent of the order in which 

the limits are taken. This is often true in asymptotic methods. 

The Aco in the first equation (4.22) refers to the bandwidth associated with the 

individual structural modes. This was implicitly true in (4.23) also, but in (4.23) a very 
large number of structural modes had been already assumed in the frequency interval of 
interest. Perhaps more significant and surprising is that in (4.22), compared to (4.23) there 
is an additional multiplicative factor, 

(fcfx.y )) 2 

(% 2 (x,y)) 

Moreover, this ratio asymptotically approaches zero as (Oc approaches infinity. This is true 
because, < X P C 2 > equals 1/2 for high center frequencies or high mode number, while <'P C >2 
becomes smaller and smaller for higher center frequencies. This implies that for large 
modal number < V F C > 2 approaches zero. Therefore, it is only the lower order structural 

modes that significantly contribute to the acoustic response. 

The above ratio appears as a consequence of assuming that the flexible wall area is 
small compared to an acoustic wavelength, such that, F r (x,y,z 0 ) on the flexible wall area 
can be approximated by F r at the centroid (x 0 ,y 0 ,Zo). From (4.18), the integral over the 

flexible area is: I vn s JjAf F i< x «y>Za)'i' in {x,y) dxdy . This assumption, removes 

F r (x,y,Zo) from inside this integral leaving only TV The entire integral is squared, hence 
the term <'F c (x,y)> 2 , (Af 2 also appears). 

Consider now two limiting cases for Af. First, if the entire wall is flexible (i.e., Af 
= the entire wall area, L x Ly), then the integral, I r,m is made up of a summation of terms 
which look like: 
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cos (r x 7 t T m x 7 t) cos (r y 7 t — m y 7t) — 1 cos (r x 7t ± m x 7t) cos (r y 7t + m y 7c) 1 

(r x 7C T m x 7t) (r y 7T - m y 7 t) (r x 7i ± m x rc) (r y rc + m y rc) 

L x Ly Lx Ly 

Assuming r » m, since there are an infinite number of acoustic modes, r, for each 

structural mode, m, the above expression is approximately: 

^ cos (r x 7t) cos (r y 7t) - 1 
(r x 7 t ) (r y 7t ) 

L x Ly 

Obviously, primarily the lower acoustical modes contribute. Therefore, the case where the 
entire wall area is flexible is not of interest here, since in this analysis a large number of 
acoustical modes are being studied. The case of the entire wall vibrating is effectively the 
case where the scale of the wall is large compared to an acoustic wavelength. For small r 
and small m, CMA should be used, of course. 

Now consider the other limiting case where the length scale of Af is much smaller 
than an acoustic wavelength. The acoustic modal function, F r (x,y,z 0 ) can be 
approximated by F r (x 0 ,y 0 ,z 0 ), yielding: 

Iyn = F/Xo.ycZo) JJ Af^ / m( x >y) dxdy. For Tm = sin(m x 7tx/L x ) sin(m y 7ty/L y ), 
the integral over T'm equals: 

L x 

(m x 7t ) 

for say, x ranging from x 0 - e to x 0 + e, and y ranging from y 0 - e to y 0 + e. Thus 
primarily, the lower structural modes will contribute in this case. The acoustic modes will 
contribute, however, for all r. Therefore, only this case need be considered for a large 
number of acoustical modes, but a finite number of structural modes, i.e. the case where 
the length scale of Af is smaller than an acoustic wavelength. 

The comparison of equations (4.22) and (4.23) shows clearly that the double limit, 
AM -» oo and AN a -4 °o, is not unique, that is, the result depends on the order in which 

AM and AN A approach infinity. In equation (22), AN a -4 «> and then AM -4 <*>. On the 
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other hand, in equation (4.23), AM->°o and then AN A This suggests incidentally, 
that the phrase Asymptotic Modal Analysis is particularly appropriate. All of these results 
using AMA are asymptotic and do break down if the assumptions under which they are 
derived are violated. It is clear, of course, that the real value of equation (4.22) is to 
demonstrate that the acoustic response of higher frequencies decreases rapidly for a large 
but finite number of structural modes, that is less in number than the acoustic modes. 

In practice this means that one must first estimate AM and AN A . The larger of the 
two is then allowed to approach infinity first and then the second. This determines whether 
equation (4.22) or (4.23) is applicable. Presumably there is yet another limit where AM = 
AN a and both go to infinity. Properly speaking, this is probably a matter of length scales, 
i.e. the dimension of the plate, the dimension of the cavity and the wavelengths of each. 
The time or frequency scales are clearly also important. 

Conclusion 

An asymptotic modal analysis result for the mean-square pressure in a bandwidth 
has been derived for the case where there are a large number of acoustic modes responding 
for each structural mode. The limit is then taken as the number of structural modes 
approaches infinity. Upon comparison of this result to previous work, where the limits 
were tak en in the opposite order, it is shown that the order in which the number of modes 
goes to infinity is important. The significant difference in the two results is a ratio made up 
of spatially averaged structural modal functions, which only appears in the newly derived 
result. This ratio approaches zero as the center frequency approaches infinity. 

Since the two approaches yield different results, it is important to determine which 
modes in the coupled system are approaching infinity faster. This information will result in 
the proper choice of an asymptotic modal analysis result for the mean-square pressure. 
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A ppendix 


In this appendix an alternative AMA expression for mean-square pressure is 
derived. Here, the variable describing the input excitation is the power spectra of point 
forces applied to the vibrating wall, rather than the power spectra of the wall acceleration in 
generalized coordinates. 

From equation (4.4) in the text and power-spectra relations: 


$c£(j(a>) = 



^m(x ,y Wx* (co;x,y,x*,y *)dxdydx*dy * 


We can relate the generalized coordinates (q's) to the generalized forces (Q's) by: 


where the transfer function Hm is given by: 




1 


M, 


2 2 
-CD +2iC m CD m O>fCO m 


It follows that: 


°^n (co) =c ° 2h ™ (co) WH„(-<D)OQ mQn (CD) 


Using the above relations in equation (4.16) in the main text: 


Ow (co;x,y,x*,y*) = X S ^ m (x,y)'I / n(x*,y*)cD^I m (cD)co^Hn(-®) 



m n 

v J / m(x,y) v f' n (x*,y*) • O^(co;x,y,x*,y*) dxdydx*dy* 


(Al) 


Equation (Al) relates the external pressure on the wall to the wall acceleration. 

The above expression can be simplified by ignoring off-diagonal coupling ("joint 
acceptance"). This means omitting all terms in the double sum except those for which 
m=n. This step is justified because Hni(CD)Hn(-CD)« IHm(co)l^ and Hn(co)Hm(-co)« 
IH n (co)|2. Additionally, d>p e is slowly varying with respect to the rapidly varying transfer 
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functions (H m (co)) and therefore can be treated as a constant independent of co. This is a 
consequence of assuming that the forcing is white noise. Equation (Al) then becomes: 

Ow (w;x,y,x*,y*) = 2 vp m(x,y)'I / m(x*,y*)co 4 |H m (co)| 2 

(A2) 


■ft ft 


X l / m(x,y) v l / m(x*,y*)* <I ) if(® m ;x,y,x*,y*) dxdydx*dy* 


The load, p e , can be described by many point forces, Fi. This is helpful in 
obtaining and understanding the results. It does not lead to any lack of generality in that 
any continuous spatial distribution may be represented by many point loads (I is the 
number of point loads). 

I 

p(x,y,t) = XSCx-x^Cy-y^FiCt) 
i = l 

The associated power spectrum is 

Cy (co;x,y,x*,y*) = 22 5 ( x ‘ x 0%-yi) 5 ( x *‘ x j*) 5 (y**yi*^ F i/ co ^ 
i j 

where O f is the cross-power spectrum of the point loads. 

ij 

Therefore, equation (A2) becomes: 

0* (co;x,y,x*,y*) = 2 v F m ( x .y)^m( x *.y*)“ 4 | H m(o>)| 2 
m 

2 2 V J / m( x i>yi) X ^ / m( x j*,yj*) 

» j 

Using this expression for O* in equation (4.15*) of the text, yields: 

Of (co;x,y,x*,y*) = II F r (x,y,z„)F r (x*,y*,Zo) 

•2 'F m (x,y) v F m (x*,y*)co 4 |H m (co)| 2 2 2 'Pm(x i ,yi) x F m (x j *,y j *) 0 Fij (co m )dxdydx*dy* 

m i i 


i J 
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This relates the power spectra of the forcing to the power spectra of the pressure. In order 
to obtain the mean-square pressure, this expression must be integrated with respect to 


frequency, by definition: 


(p<£of 


);x,y,z) = ^ 


—2 

f 2 = — = R((t=0;x,y,z) = I 0>f(co;x,y,z) dco 


Therefore, the mean-square pressure expression is: 


P 2 


(p°c?F 


1 


j | H^co) | 2 coi H m (co) | 2 dco ^ 


FXx,y,z 0 )-S Wll v Pm(x 1I yi)'F m (x J ,y J ) ^(cojdxdy ] ( A4 > 


1 J 


1 


Here, the * superscript has been dropped since the point pairs xj,yi and xj ,yj 
obey xj= xj* and yj = yj* , i.e. they range over the same points on the structure. 
Therefore, the * superscript is superfluous. Note also that the previous equation retains a 
fully modal character of the classical type. We have not yet passed from the domain of 
modal analysis to AMA. 

Figure Al, shows the smoothness of the functions IHm(co)l 2 and relative 

to IH r (co)l 2 , for the case where there are a large number of acoustic modes for each 
structural mode. The figure shows the response at one structural mode whose resonance 
frequency is co m . Implicit in this drawing is the assumption that the acoustic damping is 
less than the structural damping. In this case, IH m (co)l 2 and Op-j (co) can be treated as 

constants relative to the rapidly varying iHr(co)l 2 . 
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Figure Al. Qualitative sketch of transfer functions, IHm(co)!^ , IH r (co)l^ and power 
spectra, Op., (co). 


The frequency integral then becomes: 


HmCCOj 


2 



H/co) 


2 4 

CD dCD 


where 


HJfflJ 

2 

i 


f 2*\ 





1 


4M„C (D 


4 * 
m 


The frequency integral is equal to 


I | H/co) | 2 co 4 do) = for sma11 damping ratio £ r - 

Jo 

The higher order damping terms are not retained here. Substituting this into (A4) yields: 

iZi 


P 2 


(poC?f 


= ^ 1 ^ 1 FJMA? 


r co m 4 


V 2 


Fj(x,y,2 0 ) T m (x,y)X X ^m(xi>yi)'J'm(x>yj) dxdyj 
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This is the mean square pressure as a function of the external forcing Fjj for one 
structural mode. This expression is still of the CMA type, since there have been no 
restrictions on the summation overr. 

Assuming the flexible area is small compared to the cavity dimensions, F r (x,y,zo) 

can be approximated by its value at the centroid of the flexible area (xo.yoZo)- In the AMA 

limit, the parameters, Cr A > ®r A and M r A can be replaced by their values at the center 

frequency. The summation over r, as r — » «» can be approximated by: 

X (F I (x,y,z)) 2 »{F c (x,y,z)) 2 AN A 
r 

Therefore, the expression for the mean- square pressure in terms of the power spectra of the 


external forcing is: 

Pl_ . I CcV) 1 j 1 J 

. * A, 4 ' V 2 


(poC?f 


Mm^m C0 : 

Jj Y m (x,y )X X ^m(xiyD^m(x>yj) dxdy 


T 


(A5) 


i j 


where <Zc 2 > 0 is the ratio of <F C 2 > over the volume to <F C 2 > o ver Af . 

Equation (A5) represents the contribution from one structural mode. Implicit in this 
result is the assumption that infinitely many acoustic modes are present in the bandwidth of 
a single structural mode. The total mean-square pressure is found by summing over the 
finite number of structural modes (m) in the bandwidth of interest. 
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